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A'ba  tract 

This  report  is  the  second  chapter  of  a  proposed  'book  on 
the  Techniques  of  Applied  Matheaatics.  Because  of  the  time  lag 
■before  publication  and  because  of  the  interest  in  the  subject 
matter,  it  was  thought  advisable  to  issue  selected  chapters  as 
project  reports.  The  first  chapter  reviews  well-known  facts 
concerning  linear  algebra  and  will  not  be  issued  as  a  report. 

This  second  chapter  discusses  the  theory  of  the  Dirac 
delta  function  and  singular  functions  which  play  such  an  im- 
portant part  in  classical  applied  msthematics  and  in  modem 
physics,  especially  quantum  theory.  Although  the  results  are 
in  most  cases  well-'oiown,  it  has  been  thought  worthwhile  to 
present  p   fairly  rlgnrous  mathematical  treatment  of  these  topics 
in  order  to  give  the  physicists  confidence  in  the  use  of  the 
delta  function.  The  treatment  is  based  on  the  theory  of  dis- 
tributions which  W8  3  developed  by  laurent  Schwartz, 

Comments  and  sugrestions  by  readers  of  this  report  will 
be  most  welcome  to  the  author. 
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Introd";ction 

A  ffood  pert  of  ais themstical  physics  is  devoted  to  a  discasaion  of  the  sin- 
gtilsrlties  In  a  field;  for  ex;=mple,  point  sources,  line  sources,  dipoles ,  electric 
and  mpgnetic  currents,  etc.  All  these  phenomena  represent  soae  kind  of  discontinuity 
in  the  field  and  it  is  difficult  to  give  an  adequate  iiathemstical  description  of 

them. 

The  physicists   have   tried  to  simplify  the  npthematics   of  such  -ohenomena 

"by  introduc iiit   the  Dirac    6-function.  Its  us\ial  definition  is   this:    the   6-fui;ction   is 

zero   for  every  value  of  x  except  the  origin.,   where   it   is   infinite   in  such  a  v/ay  that 


/ 


6(x)dx  =   1. 

-00 


Mathematically,  this  definition  is  nonsense.   If  a  function  is  zero  every- 
where except  ft  one  point,  its  integral,  be  it  Riemani;,  Icheegue,  or  anyone  else's, 
is  necessarily  equal  to  zero.   fTowever,  the  6-function  has  proved  itself  to  he  of  such 
grept  utility  tha-  the  physicist  rightly  refuses  to  give  it  up.   It  is  the  mathen^a- 
tician's  task  to  find  a  method  hy  which  the  use  of  the  6-function  can  "be  justified. 

One  such  methort  has  been  known  for  a  long  time.   It  is  to  consider  the  6-functicr 
r-s  a  disgiiised  Stielt.ies  integral.   This  idea,  although  correct,  has  not  proved  very 
fruitful.   Recently,  Ifi.urent  Schwartz*  has  propose.:  a  methoa  which  not  only  justi- 
fies the  us<=  of  the  6-function  hut  also  the  use  of  all  derivatives  of  the  6-fijLnction. 
By  this  method  every  continuous  function  can  "be  differentiated  arbitrarily  often.  All 
the  stop  signals  that  the  me.thersaticians  have  imiiosed  upon  the  free  hs.ndling  of  3a,the- 
matics,  such  as  questions  about  convergence  and  the  interchange  of  li'iting  operations, 
become  much  easier  to  avoid  and  to  neglect.   The  method,  called  the  Theory  of  oistri- 
butions  by  its  discoverer,  will  be  developed  in  this  chapter. 

Interpretation  of  the  6-function 

What  should  be  understood  by  the  6-furiCtion7   Its  most  iaportrnt  property 
and  the  one  that  makes  it  so  useful  is  the  following: 

00 

(2.1)  f       £(x)^(x)dx  =  ^(0)     ; 

^-00 

this  property  may  also  be  expressed  thus:  c-iven  a  continuou<5  functio.^i  ^(x),  the  £- 

function  picks  out  the  value  of  that  function  at  thp  origin.   We  may  say  then  thrt 

^   "Theorie  des  Distributions  et  Tranformation  de  Fourier,"  Annales  de  I'Universite  de 
Grenoble,  tome  XXI,  pp.  57-7A,  19-^5;  and  "Generalisation  de  la  Notion  de  Fonction  ... 
de  Fourier,"  tome  'Lilll,   pp.  7-2^,  194.7-3.  See  also  L.  Schwartz,  Theorie  des  Distri- 
butions .  Actualites  Scientifiques  et  In^ustrielles  IIo.  1091  and  1122,  Pa^is,  1950-1. 
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the  6-f-'inctioii  Is  really  a  functional  on  the  space  of  continuous  functions  because 
it  assigns  a  number  to  every  contiaous  function. 

This  vie\'*point  may  appear  more  plausible  if  we  consider  that  every  inte^rable 
function  f{x)   may  be  interpreted  as  a.  linear  functional.   Let  G-  be  the  linear  vector 
space  of  continuous  functions  0(x)  which  vanish  outside  a  finite  interval.   The  f-func- 
tional  fissions  to  any  ftinction  j5(x)  in  the  s-osce  C  the  v?lue 

(2.2)  f(^)  =/^  f(x))«U)dx. 

This  equation  shows   thot  any  integrable  function  f (x)  defines  a   linear  functional   on 
^    .      Conversely,    sup: ose   there   is  given  a  linear  functional  on    C  ;    that   is,    to  any 
(6(x)    in  C^let   there  be  assigned  a  niimber  denoted  by  f(0)   such  that 

(2.3)  f(a0i  +   Pfig)   =  af(li^)  +   3f(02), 

where  a  and  P  are  arbitrary  numbers.   In  general,  there  will  not  exist  a  function  f (x) 
such  that  (2.2)  is  valid.   However,  if  f (0)  is  also  a  continuous  functional,  that  is, 
If  the  values  f(j!l  )  converge  to  zero  as  the  maximum  of  0  (x)  converges  to  zero,  then 
as  r,  Biesz  has  proved,  there  exists  a  function  m(x)  such  that  f (^)  can  be  written  as 
a  Stieltjes  integral: 
{2.U)  f(0)  =/"  0(x)dm(x). 

-  00 

In  case  m(x)  has  a  deriva.tive  f (x)  and  is  equal  to  the  integral  of  its  derivative, 
equation  (2.^)  becomes 

00 

f(0)  =/^   0(x)f(x)dx. 

ITow,  let  us  consider  the  6-functlonal:   To  every  j4(x)  in  C    it  assigns  the 

number  0(0).   This  definition  clearly  satisfies  (2.3)  and  thus  shows  that  the  6- 

functional  is  linear.   The  5-functlonal  also  is  continuous,  because  if  the  msjclm-^jn 

of  0  (x)  converges  to  zero,  then  surely  the  values  0  (O)  converge  to  zero.   By  the 
n  n 

theorem  of  F.  Eiesz  there  must  exist  a.  function  m(x)  such  that 

6(0)  =  0(0)  =/"^  0(^)dm(x)  . 

This  function,  m(x),  is  actually  the  Eeavlside  unit  function,  H(x),  defined  as  follows: 
H(x)  =0,       X  <  0 
=1,       X  >  0  ; 
therefore 
(2.5)         6(0)  =  /^  0(x)dH(x). 

Cf  course,  equation  (2,5.)  and  the  interpretation  of  the  6«rfunction  as  a  Stieltjes 
integral  were  well  known.  The  contribution  of  Laurent  Schwartz  consists  in  the 
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Introduction  of  the  concept  of  functionala  and  in  the  investigation  of  the  impli- 
cations of  that  concept. 

Instead  of  continuing  the  disc\ission  off  functionals  on  the  space  (\,   we  shall 
introduce  a  space  /(/  ,  which  will  he  more  useful  in  applications.   Let  ^he  the  linear 
vector  space  of  fimctions  j^(x)  which  have  continuous  derivatives  of  every  order  and 
which  vanish  identically  outside-  a  finite  interval.  Any  such  function  ^(x)  will  "be 
called  a  testing  function.   The  use  of  the  space  /(/will  enable  us  to  differentiate 
functionals  as  often  as  we  please. 

Let  T((j)  denote  the  m-unher  a  functional  T  assigns  to  any  testinr  function  ^(z). 
A  functional  is  continuous  if  the  numbers  T(0  )  converge  to  zero  whenever  the  sequence 
of  testing  funct:ons  0  (x)  converges  to  zero  in  the  sense  of  Xj\     A  seq-^ence  of  testing 

functions  0  (x)  cnnverges  to  zero  In  the  sense  of  Vif  the  mnximum  of  —  —  for  p  fixed 

d^ 
(p  =  0,1,2,...)  converges  to  zero  as  n  becomes  infinite  and  if  all  the  functions  0  (x) 

are  zero  outside  the  same  finite  interval.   It  is  an  immediate  consequence  of  this 

definition  that  if  t     converge  to  0  in  the  sense  of  ■(7  (i.e.,  ft   -0  converge  to  zero  in 

the  sense  ofc'),  then  the  numbers  T(0  )  converge  to  T(^).  Any  linear  continuous 

functional,  T,  on  the  space  </  will  be  called  a  distribution,. 

The  set  of  all  distributions  forms  a  linear  vector  space  which  we  denote  by 

^'.  This  space  is  said  to  be  the  space  dual  to  -(7  .  Any  integrable  function  f(x) 

defines  a  distribution  by  the  formula: 

(2.6)         f(0)  =/f(x)^(z)di. 

(Note:  all  integrals  will  be  from  -oo  to  +oo,  unless  otherwise  specified.)   In  such  a 
case,  we  say  the  distribution  is  a  function.   Other  examples  of  distributions  are  the 
5-functional,  which  assigns  the  value  0(O)  to  every  function  0(x),  and  the  functionals 
which  assign  the  value  0  -^  (O)  to  0(x). 

The  concept  of  distribution  is  a  generalization  of  the  concept  of  function. 
A  function  Is  defined  by  the  values  it  assumes  for  given  values  of  x.  A  distribution 
does  not  have  any  value  for  a  given  value,  of  x.   It  is  defined  only  by  the  values  it 
assigns  to  the  set  of  testing  functions.   Nevertheless,  oecause  the  function  notation 
is  so  convenient  to  use,  we  associate  a  symbolic  function,  t(x),  with  every  distribution 
T,  and  write 

(Z.-^)         T(/5)  =/t(x)0(x)dz. 

The  function  t(x)  and  the  Integral  in  (2.'')  will  have  no  meaning  except  that  given  by 
the  left  side  o'  (2.7).  However,  if  T  is  defined  by  an  integrpble  function  f(x),as  in 
(2.6),  we  shall  identify  t(x)  with  f(x). 
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With  the  6-functional  we  aBSociate  the  6-f-aiictlon,  6(x),  which  is  p  sym- 
holic  function  and  has  no  values  or  meaning  othtr  than  what  is  implied  "by  the  formula 

0(0)  =/6(x)0(x)di. 

Note  that  since  a  distribution  is  8.  continuous  functional, 

lim  T(f!S  )  =  T(lim0  ). 
n         n  ' 

where  the  lim  0  is  to  he  understood  in  the  sense  of  convergence  in  </  ,   In  terms  of 
n 

symholic  functions  this  equstion  can  he  expressed  as  follows: 

Theorem  I.   If  t(x)  is  a  symbolic  funct-Jon  and  t>    i^)   is  a.  sequence  of  testing  functions 


converging  in  >v  ,  then 

'lim/t(x)^  (x)dx=/t(x)  lim  0  (x)dx. 
Prohlems 

2.1  Prove  the  statement  in  the  text  that  an  integrahle  function  f (x)  defines 

a  distribution  hy  the  formula  (2.6).   Hint.   Show  that  /f (x)0(x)dx  is  a  linear  contin- 
uous functional  on  /(/  . 

2.2  Show  that  the  function 

2^      1 

0(x)  =  e'^^e"^^-^)''  ,0<x<a 

=  0  ,x<0,ora<x 

is  a  testing  function. 
Translations  of  Distrihutiona, 

We  shall  now  develop  a  calculus  of  distributions  and  their  associsted  sym- 
bolic functions.   Our  aim  is  to  extend  to  distributions  the  ideas  and  concepts  which 
are  valid  for  ordinary  functions.  This  extension  is  usually  obtained  by  re-writing 
these  ideas  and  concepts  so  thpt  they  apply  to  functionals. 

Consider  the  concept  of  the  translation  of  a  distribution.  Supi  ose  the  dis- 
tribution is  defined  by  a  function  f (x);  then  the  translation  of  it  through  a  distance 
h  is  f(x+h).  and  the  corresponding  functional  is 

(2.8)  /f(x+h)^(x)dx  =/f(i)0(x-h)dx. 

Write  fyf   for  the  translation  of  f(x)  thro-igh  a  distance  h:  then  (2.8)  becomes 

(2.9)  t^f(t)   =  f(t^i)   . 

Equ-tlon  (2.9)  is  used  to  define  the  translation  of  a  distribution.   Let  T  be  an  arbit- 
rsry  distribution  and  let  "t  T  'be  its  translation  thro-ugh  an  amount  h:  then  we  define 

(2.10)      'i^^^iii)  =  -r{'i^_J>). 

The  corresponding  symbolic  function  will  be  written  t(i+h),  and  from  (2.10)  we  have 
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/t(3t+h)^(x)dx  =/  t(i)0(x-h)dJt:; 

thus   t(x)  actt\a.lly  liehaves   like  a  fxinction. 

Suppose  T  is  the  6-functlonal;  then 

/6(i+h)0U)di  =/ 5(x)0U-h)djc  =  0(-h); 

consequently  d(rfh)  Is  the  symholic  function  for  the  distribution  which  assigns  the 

value  0(-h)  to  every  testing  function  0(x). 

Differentiation. 

The  derivative  of  a  distribution  T  is  defined  in  a  manner  similar  to  the 
definition  for  a  function: 

r,T(0)-T(0) 

(2.11)        T'(0)  =  lim  -^ . 

h-^ 

From  (,?.10)  we  get 


T«(0)  =   lim[T(r  j^0)  -  T(0i]/h  =  lim  t[^^=| J 

=   nm  T    [iC-Hlz^    =  T    [lim  ^i-^h^i-)^ 


=  T(-0«(x))   =  -T(0»)    , 

where  use  has  been  made  of  the  fact  that  T  is  a  continuous  functional.  We  state  this 

result  as 

Theorem  II.   If  T  is  a  di3tri"hution  and  j}  is  a  testing  function,  then 

(?.12)        T'(0)  =  -T(0')  . 

let  t(x)  be  the  symbolic  fxmction  associated  with  T  and  let  t'(x)  be  the  symbolic  func- 
tion associated  with  T'jthen  (2.12)  takes  the  form 
(2.13)        /t'(x)0(x)di  =  -/t(x)^'(x)dx, 

which  is  exactly  the  result  that  would  be  obtained  if  t(x)  were  a  real  fxmction  and 
integration  by  parts  were  used. 

To  illustrate  the  concept  of  differentiation,  consider  the  Heaviside  unit 
fimctlon  H(x).  Since  it  is  a  summable  function,  it  may  be  considered  as  defining  a 
distribution.  As  such  what  is  its  derivative?  By  (2.12) 

H«(0)  =  -H(0')  =  -/^   0'(x)dx  =  0(0)  , 

since  0(x)  vanishes  for  sufficiently  large  x;  therefore,  we  may  write 

(2.1i*)  H'  =  6 

since  H'  assigns  the  value  «5(0)  to  every  testing  function  0(x).   Note  that  the  deri- 
vative of  H(x)  considered  as  a  function  is  zero  everywhere  except  at  the  origin,  where 
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it  is  undefined,  while  the  derivative  of  the  distribution  defined  by  H(x)  is  the 
6-f'anctional. 

What  is  the  meaning  of  6 '7  3y  (2.12) 

(2.15)  6'(«5)  =  -6(0')  =  -|*'(0); 

hence  the  6 '-functional  assigns  the  value  -0*(O)  to  every  testing  fimction  0(j').   In 
the  same  way,  we  can  prove  that 

(2.16)  6(^\0)  =  (O'^^^^^O). 

Sun-pose  T   is   defined  by  en  absolutely  continuous   function  f  (x)^  thfi.t   is,  a 
function:  which  is   the   intp^^'ral  of  its   derivative;   then  by  integration  by  parts 

-7'f(x)0'(x)dx  =yf '(x)0(3:)d2j 

consequently  we  find  that  the  symbolic  function  for  the  derivative  of  the  distribution 
defined  by  f (x)  is  the  spjne  ps  the  derivative  of  f (x)  considered  as  a  function.   However, 
if  f(x)  is  not  absolutely  continuous,  the  two  derivatives  sre  different.   Henceforth, 
the  derivative  of  f(x),  written  f '(x),  will  mean  the  symbolic  function  for  the  deriva- 
tive of  the  distribution  defined  by  f(x). 

We  shall  prove 
Theorem  III.   Suppose  f (x)  is  absolutely  continuous  and  differentiable  except  at  the 
points  a  .  a^. ..  ,a  .  where  f(x)  has  .l^amps  of  magnitude  g..,  a^  ....a^^  respectively;  the. 

f'(x)  =  [f '(x^l  +  a^6(x-a^)  +  a.^h(,7.-B.^   +  ...  +  a^6.(x-a^) 
whfTe[f '(xJ3  is  the  derivative  of  f (x)  considered  as  a  function. 
To  prove  this  theorem,  put 

(2.17)  g(x)  =  f(x)-a^  H(x-a^)  -  a2H(x-a2)  -  ...  -a^H(x-a^). 

It  is  clear  that  e(x)  will  be  absolutely  continuous  and  differentiable,  so  that  the 
derivative  of  g(x)  will  be  the  same  whether  g(x)  is  considered  a  distribution  or  a 
function.   When  (2.1?)  is  differentiated  as  an  equation  involving  distributions,  we 
find  that 

(2.18)  f '(x)  =  g'(x)  +  a^6(x-a^)  +  CL^kiji-&^   +  ...  +  a^£(x-a^). 

Since  g'(x)  is  the  same  as  the  derivF.tive  of  f(x)  considered  as  a  function,  we  see 
that  the  derivative  of  a  function  considered  as  defining  a  distribution  is  the  us-ual 
derivative  t)1us  the  sum  of  6-functions  at  the  jumps,  multiplied  by  the  magnitude  of  the 
jumps. 

In  a  similar  way  we  see  that 
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(2.19)    f"M  =   P"U)  +  a^6'(x-a^)  +  ...  +  a  6  •  (x-a  ) 

+  e^6(x-b^)  +  ...  +  Pj^6(i-bj^), 

where  b,  ,  'b^,'*»,'b,    are  the  points  of  discontinuity  in  g'(x)  and  P,,  3  ,.«•,?,  sre 
the  corresponding  .iumps  in  the  \-alue  of  g'(x), 

Examples  of  Derivrtives  of  Distrihutions 

Consider  the  distribution  defined  by  the  function  |x|.   Its  derivative  is  the 
function  which  is  -1  for  the  x  negative  and  +1  for  x  positive.  We  call  this  function 
signum   x  and  denote  it  by  sgn  x.  We  have 

(2.20)      d         _.,  V 

—  sgn  X  =  26 (x)  , 

since  the  sgn  function  is  constant  except  for  a  jump  of  magnitude  2  at  the  origin. 

As  an  application  of  this  result,  consider 
.1 
I  =y    |x|v)/«(x)dx  . 

Using  repeated  integration  by  parts,  we  find  that 

I  =  lxU'(x)I_^  -/_^  lx|«^4,'(x)dx 

I     1     1 

=  [Ix!  ^V^M    -    Ul  v;,(x)]_^  +/^  |x|%(x)dX 

=  u/«(i)  -  ./(-I)  -  y,(i)  -,^(-1)  +  aj>(o)  . 

This  method  is  simpler  than  the  usual  one  where  I  is  written  as 

■y  aD(»"  (x)dx  +  y  x)/"(x)dx 
1  0 

and  integration  by  parts  is  applied  to  both  intpfTsls. 

As  another  illustration,  consider  the  problem  of  maximizing  the  integral 

K  =y_^  f  (x)(x-Xj^)(x-i2)dx, 

where  f (x)  is  an  integrable  repl  function  such  that  |f (x)|  <  1,  and  where  x^,  x  are 
arbitrary  numbers  between  -1  and  1.  For  given  x^,  x  ,  K  will  be  a  maximum  if  f (x)  =  +1 
when  (x-x^)(x-X2)  is  positive  and  if  f(x)  =  -1  when  (x-x^)(i-X2)  is  negative;  therefore 
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K<K*=y_^  (x-x^)(x-i2)sgn[j;x-x^)(x-X2D'i^  • 

Now  K*  wil]  "he  a  msxlmtm  for  those  values  of  x^   and  x^   which  make 

We  have 

||i  =  .y_^U-x2)8rnn(x-x^)(x-X2[ld^ 

1 
+/|^  (x-x^)(x-X2)26(x-x^)dx  . 

The  second  inte#:ral  is  zero,  since  (x-x,)(x-x_)  vanishes  when  x=x^.  finally, 


9g* 

ax. 


^   (x-X2)dx  +  A   (3!:-X2)di  ~J^   {x-x^)&i. 


ax* 

There  is  a  similar  formula  for  r —  .  The  remainder  of  the  solution  niay  be  left  to 

axg 

the  reader. 
Pro> 1 ems 

2.3  If  f(x)  is  absolutely  continuous,  show  that  the  derivative  of  the  dis- 

tribution defined  by  f(x)H(x)  is  f«(x)F(x)  +  f(o)6(x). 

2.^  If  f(x)  is  sbsolntely  continuous  with  zeroes  at  x, ,...,x  ,  the  deri- 

i      n 

vative  of  sgn(f(x))  is 

26(x-x^)3gnf' (x^)  +  ...  +  26(x-x^)sgnf '(i^). 

2.5  Shov*  that  the  integral 

u(x)  =/  lx-tl(«(?)dt 
satisfies  the  differential  equation 

u"  =  20 (x). 
Hin*-..  Use  (2,20) 

2.6  Show  that  the  integral 

u(x)=/e^^'^-^l0(OcLC 
satisfies  the  equation 

u"+  ic^u  =  20  (x). 

2 
Hint.   Use    (2.20)  and  the  fact  thnt    (sgnx)     =  1. 

2.7  If  p(x)    is  an  integrsTjle  non-negative   function,   show   that 

is  a  solution  of 

(pu-)'  =  £(x). 
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F^-nctions  with  Infinite  Jumps 

What  is  the  derivative  of  the  distribution  corresponding  to  the  function 

-1/2 
which  is  zero  for  negative  x  and  x  '   for  positive  zT  The  usual  derivative  is  zero 

— '^/? 
for  nep- tive  x   and  -x"'  /2  for  positive  x.   However,  the  letter  function  is  not  in- 

tegrahle  and  cpnnot  define  a  distribution. 

\'Ie   should  like  to  use  formula  (2,12)  to  define  the  derivative,  hut  integra- 
tion by  riarts  cannot  be  applied  to  the  integral 

Vu)i"^^^<ix  . 


-/ 


The  solution  is  to  realize  that  this  integral  is  the  limit  as  c  — >  0  of  the  following 
integral: 


-^ 


0«(x)i"^/2dJC. 


Now  integration  by  parts  is  possible  and  we  have 

-/^'(x)x"^/^dx  =  -Um   /^  0'(x)x"^^^dx 
t-^  0 

(2.21)  ,      t»      , 

=  lim  [iS(c)c-^/2  l/'^(^)^-3/2^-|  ^ 

€->0  '^ 

This  last  expression  was  first  considered  by  R'^damp.rd.   He  called  it  the  finite  part 
of  the  Integral 

We  generalize  this  concept  to  a  function  g(x)  which  is  not  integrrble  because 


of  a  di?continuity  at  x=a.   Consider  the  integral 

lb 

g(x)dx  . 


/ 


This  integral  may  have  an  asymptotic  expansion  in  positive  and  negative  powers  of  c. 
The  terms  in  negative  powers  of  e  or  in  log  e  are  called  the  infinite  part  of  the  in- 
tegral.  The  limit  of  the  difference  between  the  integral  and  its  infinite  part  as  € 
approaches  zero  is  called  the  finite  part  of  the  integral.  This  finite  part  of  the 
integral  will  be  used  to  define  distributions  corresponding  to  non-integrable  f-onctions. 
An  important  set  of  distributions  is  obtained  in  this  way  when  we  take 

g(x)  =  XT  where  the  real  part  of  p  <  -1.  Assume  first  that  p  is  not  a  negative  integer 
and  consider  the  finite  part  of  the  integral 

.00 

J     xP0(x)dx. 
o 

Suppose  that  m  is  the  smallest  integer  such  thrt  the  real  part  of  (p+m)  >  -1;  then  by 

using  the  Taylor  series  for  0(x)  it  can  be  shown  that  the  following  integral  exists 
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In  the  ordinary  sense; 


(m-1), 


1^(0)  =/xP[i5(xM(o)-z^'(o)-|j-^'(o)-  ...  -f^iryr-^^^^^oDi^ 

00  m— 1  /     ,  \ 

(2.22)        =  lim    ;^^I]5(xM(0)-xi!5'(0)-   ...   "  f^I^Jj"     ^"^(02    dx 


€    ->0 


=  lim 
c  ->0 


/    xJ?0(x)dx  +  '"yi         +        V2         -^     21       ;^-  ••■'    (m-lji        Ti^ 


This  shows  that  I  (0)  is  the  finite  part  of  the  integral 

/°°3?0(x)dx. 
0 

If  X  pees  from  -00  to  0,  let  us  consider  the  finite  part  of  the  integral 
y  xP0(x)dx. 


Tt  will  "be 

(2.23)        I_(0)  =/Vr0(x)-0(O)-  ...  -  ^^  0^"^-^^(O)]  dx. 

-00  ^ 
We  shsll  use  the  symhol  ('x)^  to  represent  the  symholic  function  for  the  dis- 
tributions whose  values  are  given  hy   I_|_(0)+I_(0);  that  is, 

(2.2i*)        /(•x)P0(x)dx  =  I^(0)+I_(i/5). 

The  symholic  function  for  1^{^)   will  he   denoted  ^7  (•x)^H(x);  conseq.uently 

(2.25)        /(•x)PH(x)^(x)dx  =/^('x)^0(x)dx  =  1^(0). 

Similarly,  ('x)^H(-x)  will  denote  the  symholic  fiinction  corresponding  to  I_(0). 

If  p  is  a  negative  integer,  say  p  =  -m,  then  the  finite  part  of  the  integral 

/°°xP0(z)dx 


will  he 
lim 


00  -m+1        — art-2 

/^x-0(x)dx-0(o)^^::^-0(o)^^::^ 

€->0  '- 


-u:T)i     ^"^^  '] 


(2.26) 


=  lim  /V-'"[0(x)-0(O)-x0'(8)-...-f^^  ^(^-^^0)]  dx 


-  82  - 


=  /x-"[0(x)-0(o)-T^'(O)-...-|^^0('°-^^O)Jdx 


0U)-^(O)-x^'(O)-...-  ^^^  0^"-2^(o)^dx  . 


As  an  application  of  these  ideas  consider  the  derivative  of  the  dlstrihutlon 
defined  "by  the  integrpTile  function  logjil,  Using  Theorem  III,  we  have 


-/log|x|^'(z)dx  =  -  lim  F/"   +/   |log|x|0'U)dx 

=  lim  riog€(^(cM(-o)  +f/'lvr  )^  ^1 

(2.27)  c->0l-        V  J      W-oo/c    J      X  J 

=  -[/::<"]  ^-. 

hecanse  since  the  testing  function  ^(x)  is  continuous  ajad  differentialjle  at  the  origin, 

the  term  erising  from  the  integration  iDy  parts  is  zero.  The  last  expression  in  (2.27) 

is  customarily  called  the  princlTiel  value  of  the  integral  and  is  represented  by  the 
symbol 

We  shall  use  ('x)~  to  represent  the  symbolic  function  for  the  distribution  whose  values 
are  given  by  this  expression;  then 

(..as,      /^^.P.../^^=n«  [/-%/;]  ^.. 

€— >0  *-  -" 

Equation    (2.27)   now  shows    that 

(2.29)  (logUI)'  =   Cx)-^  . 

By  the  use  of  (2.29)  we  may  find  the  derivative  of  log  x.  But  first,  log  x 
must  be  defined  for  all  values  of  x  from  -  oo  to  +oo.  For  x  positive,  log  x  =  log|i|, 
while  for  i  negative  we  assume  log  x  =  log|x|+lTi.  We  may  then  write 

(2.30)  log  X  =  log|x|  +  iTT  H(-x) 
and  then 

(2.31)  (log  x)'  =  77  -  iTT6(x)  , 

a  formrila  already  given  by  Dirac. 

Problems 

2.8  If  0(x)  is  a  testing  function  prove  that 
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fii 


ProTjlem  2.9        If  a  >  0,   prove  that 


y'     / ,    vp-1  -ax  Rp) 

0    ('^>        «        -^  =        p 


aP 

for  all  values  of  p,  except  p  =  0,  -1,  -?,... 
Problem  2.10  Prove  that 

^  C'x)PH(xa  =p('x)^Vx) 

if  p  is  not  a  negative  integer. 

Prohlem  2.11.   If  m  is  an  integer  prove  that 

Problem  2.12.   If  m  is  an  integer,  prove  that 

I-  ('x)-"'  =  -anCx)-'"-^  . 
ox 

MTiltiplication  of  Distributions 

In  general,  two  arbitrary  distributions  cannot  be  multiplied.  For  example, 

o 
consider  the  distribution  associated  with  the  integrable  function  f(x).   If  f (x)   is 

not  integr^ble,  then  the  product  of  the  distribution  with  itself  is  not  a  distribution. 

One  case  in  which  we  are  certain  the  product  will  alwa/s  exist  as  a  dis- 
tribution is  when  one  of  the  factors  is  defined  by  an  analytic  function.  Let  w(x) 
be  an  analytic  function.  We  define  wT  as  follows: 

(2.?2)        vT(^)  =  TM). 

lote  that  if  vv  is  analytic  and  0  is  a  testing  function  then  \i/j!i  is  also  a  testing  func- 
tion and  so  the  right  side  of  (2.32)  always  has  a  meaning.  The  reasons  for  the  defini- 
tion (2.32)  can  be  seen  if  we  assume  T  is  defined  by  a  fijinction  f  (x);  then 

i,/f(0)  =/i/(x)f(x)0(x)dx  =/f(x)v(/(x)0(x)dx  =  fim^)    , 
which  is  the  same  as  (^.32). 

We  shall  illustrate  the  definition  by  some  examples.   Consider  x  .  —  , 

By  (2.32) 

X   .  ;i-(^)  =.^(x?i)=  lim  [/^+/*]0(x)dx  =  /0(x)dx  =  l(^); 

consequently 

As  another  example,  consider  x6(x).   By  (2.32) 
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x6(^)  =  6(i0)  =  0, 

since  the  value  of  x0(x)  at  the  origin  is  zero.  We  may  write 

(2.33)        x6(x)  =  0. 

Even  though  the  product  was  defined  only  when  one  factor  is  en  analytic  function, 
it  is  possible  in  some  cases  to  extend  the  definition.   Suppose  f (x)  is  a  function  con- 
tinuous at  the  origin;  then  we  may  define  f(x)6(x)  hy  using  (2.32): 

f6(^)  =  6(f0)  =  f(O)0(O); 
and  we  see  thet 
(2.3^)        f(x)6(x)  =  f(0)6(x). 

The  derivstive  of  a  product  \(/T  will  follow  from  (2.32)  and  from  the  definition 
of  derivrtive.  We  have 

(wT)«(^)  =  ->.;T()ft>)  =  -T(vi/0«)  =  -f(n;^) '->)/' S 

=  -T[j:,i/0)i]  +  T(v,/«0)  =  T'(vv^)  +  T(h('0) 

=  (m/T'  +  ^v'T)(j^S), 

and  thus  we   o'htain   the  usual  theorem  for  differentiating  a  product: 
Theorem  17.      If  T   is  a  distribution  and  u/  a  testing  function,    then 
Oj/T)'  =  lyT'  +  >y>T, 
If  this   is  applied  to    (2.33)  we  get 

0  =    [x6(x[l  •  =   6(x)  +  x6'(x) 
or 

x6'(x)  =  -6(x)  . 
The  n-th  derivative  of  (2.33)  gives  the  general  result: 
(2.35) 


x6('^^(x)  =-n6^'^^^(x)  . 


It  should  he  noted  that  in  certain  cases  multiplicp.tion  may  not  he  asso- 
ciative. For  example  hy  (2.33) 

7^  •  (x6)  =  0 


while 

Problems 
2.13 

2.1^ 


(  -j^)  .6=6. 

If  q  >  0  prove  that 

x^Cx)P=  Cx)^^  . 

If  n  and  m  are  positive  integers  prove  that 
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a)  x""  6^"'^x)  =  {-)\l   6(x)  ; 

b)  x'^^^^'^x)  =  0; 

Hint.  Use  (2.30).  Note  6^°^(x)  is  6(i),  S^'^'M   is  6«(x),  etc. 
Problems 

2.15  If  p(x)  is  a  function  which  is  continuously  different iable  at  the  origin, 
show  with  the  help  of  (2,30)  that 

p(x)S'(x)  =  -g«(0)6(x)  +  g(0)6'(x). 
Compare  this  result  with  that  obtained  "by  differentiating  both  aides  of  the  identity 
g(x)6(x)  =  g(0)6(x). 

2.16  Show  that  the  results  of  Problem  2.10  and  2.11  may  be  obtained  formally  by 
differential  in,?  ('x)~'"h(x)  as  a  product  and  usinp  Problem  2.1%. 

Convergence  of  Distributions 

A  sequence  of  distributions  T  is  said  to  converge  to  zero  in  the  weak 
sense  if  the  numbers  T,  (|6)  converge  to  zero  for  every  testing  function  0(x).  The 
sequence  T  converges  to  zero  in  the  strong  sense  if  the  numbers  T  (0)  converge  to 
zero  for  all  testing  functions  ^(x)  and  if  the  convergence  is  uniform  for  all  test- 
ing functions  which  are  zero  outside  the  same  fixed  interval  and  have  their  deriva- 
tives uniformly  bounded.   It  can  be  shown  that  the  two  definitions  are  equivalent*. 
We  shall  usually  v^ork  with  the  first  definition. 

The  sequence  of  distributions  T  will  be  said  to  converge  to  the  distri- 
bution T  if  the  sequence  of  distributions  T, -T  converges  to  zero.   If  we  express 
the  definition  of  the  convergence  of  distributions  in  terms  of  the  corresponding 
symbolic  fimctions,  we  have  the  following  useful  result: 

Theorem  T.   If  t.  (x)  is  a  sequence  of  symbolic  functions  which  converges  to  t(x)  and 
if  ^(x)  is  any  testing  funct lengthen 

lim/t^(x)0(x)dx  =/ t(x)^(x)dx. 

Of  course,  the  convergence  is  that  for  distributions,  so  this  theorem  is 
really  just  a  re-stntement  of  the  definition  of  convergence.   Suppo^ie,  however,  that 
the  symbolic  functions  t  (x)  were  actual  functions,  as  they  wo^ild  be  if  f'^.e  correspond- 
ing distributions  were  defined  by  Integrable  ftmctions:   How  would  the  limit  of  ty(x) 
In  the  ordinary  sense,  compare  with  the  limit  in  the  sense  of  distributions?  There  are 

•  L,  Schwartz,  Theorie  des  DistributionB  ,  op.  olt.,  tome  1,  p. 73. 
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many  poasitilities :  first,  the  limit  in  the  sense  of  ordinary  point-wiBe  conver- 
gence of  t,  (x)  nay  not  exist  while  the  distribution  limit  may  exist,  second, 
the  two  limits  may  exist  hut  "be  different;  third,  the  limits  may  e::^ist  and  "be 
the  same. 

As  an  illustration  of  the  first  poS8i"bility  consider  the  limit  as  a— >  0 
of  the  set  of  functions 

t(x,a)  = — r-  . 

TT(x^+a2) 

For  X  /  0,  the  limit  of  t(x,a)  In  the  ordinary  sense  la  zero,  while  for  x  =  0  the 

limit  in  the  ordinary  sense  does  not  exist.  We  shall  prove  that  in  the  sense  of 

dlstri'butlons, 

(2.36)        lim   1— r-  =  6(x). 

a->0  tt(x  +a^) 

let  jj(x).  "be  a  testing  function.  Consider  the  integral 


(2.?7) 


-/ifarv^^^;^-^^/ 


aox 


TTdSa'^)   "   TiCx^+a-^)         ^   '  x^+a^ 

r  /-I   .00   al3(x)-|Z5(0i]       .^a0i(xM(0i] 
=  [/   +/   ]  2~~2 ^^J  2  ? dx+0(o) 

for  any  value  of  >^  , 

Since  ^(x)  is  continuous  at  x  =  0,  there  exists  a  num'ber  >^  such  that 

|0(x)  -  0(0)1  <€ 
for  any  ar"hitrary  value  of  c  >  0  and  for  -\<  ^  <.\\   therefore,  the  second  inte- 
rral  on  the  ri?-ht  of  (2.??)  can  "be  estimated  as  follows: 

,jal0(x)-^(O)l      .ac  /-l   dx   .a.  /~    dx 


«/-*   /  2^  2^        -  TT  «/ „   2  2    -n  J - 
\   TTlx  +a  )  'L  X  +a 


00  2_^  2 
X  +a 


Since  ^(x)  is  a  testing  function,  it  is  bounaed  for  all  x.  Let  M  "be  a  "bound  for 
0(x);  then 

|0(x)  -  ?5(0)1  <  2M 
for  all  X  and  the  first  integral  on  the  right  of  ^Z.V)   is  less  In  absolute  val\ie 


> 

than 


00 


-00       tt(x  +a  ; 


—  arccot  -fc  =  —T  arctan  ~ 
17        a   -rf        X 


Once  the  value  of  ">1  is  determined,  a  can  be  chosen  small  enough  so  that 


—  arc  tan  =•  <  c  . 
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Consequently, 

|M(0)!   <  2c 

for  all  values   of  a  which  are  sufficiently  near   to   zero.      Since  e  can  te  made 
arbitrarily  smell,    it  follws   that 

lim     I=lim     f'-^^       =H0)   . 
a->0  €->0''    TT(x^+a    ) 

This  result  proves (2. 36). 

The  second  possilDility,  that  the  limit  in  the  ordinary  sense  may  exist 
hut  he  different  from  the  limit  in  the  sense  of  distributions,  is  illustrated  by 
the  derivptives  of  the  set  of  functions  given  above,  i.e., 

t(x,a)«  =  - — —     . 

For  all  values  of  x  the  limit  of  t(x,a)'  in  the  ordinary  sense  is  zero,  but  we 
shall  prove  thrt  the  limit  in  the  sense  of  distributions  is  6'(x). 

The  effect  of  t(x,a)'  on  any  testing  fiinction  is  fo'iond  as  follows; 

Tr(z  +a  ;       tt(x  +a  )  tt(x  +a  ) 


a[]5(x)-^(0[! 
tt(x  +a  ; 


Tr(x  +a  ; 


In  the  computation  we  have  used  the  fact  tliat 
xdx 


/ 


U^+a^)' 


=  0. 


Since  lim  1(0')  =  0'(O),  we  have 
a->0 

lim  J(0)  =  -lim  1(0')  =  -0'(O)^ 
a->0        n->0  ' 

and  consequently  we  have  proved  that 

lim   t  (x,a)'  =  5'(x). 

The  third  possibility,  that  the  limits  in  both  senses  are  the  same, 
will  occur  in  the  case  stated  in 

Theorem  71.   If  a  sequence  of  integrable  functions  t.^(x)  converges  in  the  ordinary 
sense  to  a  function  t(x),  and  if  the  sequence  t  (x)  is  bounded  by  an  inte,s"rable 
function  f(x),  l.e, 
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\t^{x)\  <  f(x),  k  =  1,  2,  ... 

then  the  limit  of  the  sequence  of  dl5trl"butlon3  defined  ty  t.  (x)  Is  the  distrroution 
defined  by  t{x). 

Consider  the  dlstrihutions  defined  by  t,  (x).  These  are  the  Integrals 

\ii!>)   =/tj^(x)f!5(x)di. 

We  now  use  a  well-lcnown  theorem  on  Lebesgue  integration*  which  states  that  under  the 
hyj^othesis  of  Theorem  VI 

lim/  tj^(x)0(x)dx  =/t(x)0(x)dx. 

This  result  implies  that 

lim  T^(^)  =  T(0) 
and  so  proves  the  theorem. 

In  the  discussion  of  the  second  possihility  we  proved  that 

lim  —  t(x,a)  =  —  lim  t(x,a). 
a— >0  a— X) 

This  result  is  no  coincidence.  We  shall  see  that  it  follows  from 

Theorem  VII.   If  a  sequence  of  distributions  T  converge  to  a  limit  distrlTiution  T, 

then  the  sequence  of  derivatives  T'  converges  to  ? ' ,  the  derivative  of  the  limit . 

3y  the  definition  of  the  derivative  (2.12), 

T^^(i«)  -  Tt(0)  =  T(0')  -  T^(0'). 

Since  T  converging  to  T  means  that  T  (0)  -  T(0)  converges  to  zero  for  any  testing 
function  (i(x),  and  since  ^'(x)  is  a  tpsting  function  whenever  (){x)   is,  it  follows 
thnt  T,'(j4)  -  T'(0)  converges  to  zero  for  any  testing  function  ?5(x).   This  implies 
that  T'  converges  to  T'  and  so  proves  the  theorem. 

Since  the  limit  of  the  derivatives  converges  to  the  derivative  of  the  limit, 
we  may  apply  Theorem  VII  to  the  derivatives  and  conclude  that  the  limit  of  the  second 
derivatives  is  the  second  derivative  of  the  limit.  Of  course  this  process  may  he  re- 
peated. We  have,  therefore,  the  following 

Theorem  VIII.   If  a  sequence  of  distrlhutions  T  converges  to  a  limit  distribution  T, 
then  the  sequence  of  n-th  derivatives  T.:  converges  to  T   ,  the  n-th  derivative  of 
the  limit. 

This  result  cpn  also  be  applied  to  convergent  series  of  distribution.  A 
series  of  distributions 

S,  ^S^*  ... 

will  >'p   snirt   tn  conv^rp-e   if  the  partial  sum 

\=^1^^2"    •••^\  

•  Titchmareh,  Theory  of  runetlons.  page  3^*5.  The  Clarendon  Frees,  Oxford,  1932 
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convf-r^ps  in  the  sense  of  distributions.   V^e  may  now  re-state  Theorem  VIII  so  that  it 
applies  to  series: 

Theorem  IX,  If  a.  series  of  distributions 
S^  +  S^  +  ... 

conTrer|g:es  in  the  sense  of  distributions  to  a  limit  T,  then  the  series  of  n-th  derivatives 

sM  .  s(»>  .  ... 

converpes  to  j^ 

As  en  illustratidin  of  these  theorems,  consider  the  function  loe(x+ia)  where 
a  is  a  nnn-negativp  number.  We  have  by  (2.30) 


lim  lof-(x+ia)  =  log|-3c|  +1ttH(-x) 
a-X) 


and  then  by  (2.31)  and  Theorem  711 

1 


lim 


=  lim 


a->0  "^^^   a->0 


2^  2 
X   +a 


-  i 


2^  2 
X  +a 


-iTT6(x). 


This  formula  proves  afain  that 


^^"^   — 2~T' 
a-X)  tt(x  +a^) 


=  6(x) 


and  proves  also  that 


a— >0  X  +a 

a  result  which  we  have  obtained  in  Problem  2.8. 

The  theorems  of  this  section  show  thrit  when  we  work  with  distributions  or 
symbolic  functions  we  ma.y  without  danger  interchange  the  operation  of  ta^^inr  the  limit 
with  the  operations  of  integration  or  differentiation.   We  now  prove  a  theorem  which 
shows  that  we  may  differentiate  under  the  integral  sign. 
Theorem  X.   If  0(x.\)  is  a  testing  function  of  y   for  a-11  values  of  \  in  a  neighborhood 


~wi:. 


of  X  .and  if  r^  is  a  testing  function  for  \  =  X  .  then 

-o-> ex = o^ 


dr  ^  («5) 


) 


x=x. 


0  '         0 

or  expressed  In  terms   o'f  3ym>'olic  functions, 


f^/t(x)ji(x,X)dx=/t(x)|^|,     dx. 


From  the  definition  of  a  derivative,  we  have 


tx'<^> 


-lim 
X=X  AX->0 

0 
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^^X 


=  lim 


r0(x,X  tA> 

4 — k- 


,tAX)-0(xAg) 


■1 . 


Now,  liecause  a  distribution  is  a  continuous  function,  or,  what  is  eq.uivalent,  'because 
of  Theorem  I,  we  aarj   take  the  limit  inside  the  distribution  and  get 


k'  (^) 


=T 


X=X 


lim 


0(x,X^+4X)-0(x,X^) 


^X 


^    ^dX^  X=X 


This  completes  the  proof  of  the  theorem. 

We  shall  apply  Theorem  X  to  f:et  a  result  which  will  he  useful  when  we  apply 
the  theory  of  distributions  to  the  solution  of  differential  equrtions.  We  start  with 
a  simple  result  from  the  theory  of  functions  of  a  complex  variable.   Consider  the 
contour  integral 


^=/s 


Ikx 


dk 


S   k 
where  3  is  the  positively  described  semi-circle  of  radius  R  in  the  upper  h;;.lf  plane. 

We  shpll  use  Theorem  71  to  obtain  the  limit  in  the  sense  of  distributions  of  L,  as  E 

ie 
goes  to  infinity,  5'irst,  suppose  x  >  0  and  put  k  =  Re   ;  t)ien 


I^  =  l/"      exp(ixEe^®)de  . 


How, 


i3> 


Ieip(ixRe      )|   <  e 
therefore  we  see  that 


-Bis in© 


tt/2      . 
1^1   ^2/         e- 


■Rxsine 


de. 


It   is  easy  to  show    (see  Problem  2.18)    that  for  0  <  ©  <  Tr/2, 
sin  9  >  ZQ/tt; 


consequently  we  have 


■n/2 


1,1  <  2/    '.  -^x«/",e  =  ^  (l-e-^-), 


and  this  shows  that  for  x  >  0  the  limit  in  the  ordinary  sense  of  I^  is  zero.   Second, 

iQ 
for  X  =  0,  again  put  k  =  Re   ;  then 


h  = 


=   ni. 
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Finally,  for  t  <  0  the  contoiir  may  "be   replaced  "by  the  semi-circle  S'  of  radius  R  in 

the  lower  half-plane  if  we  take  into  accoiint  the  residue  at  the  origin.  We  have  then 

.    ikx 
^  =  -/s«  ^k ^*   residue. 

The  integral  over  S»  approaches  zero  for  x  <  0  in  the  same  way  that  the  integral  over 
S  approaches  zero  for  x  >  0.   Since  the  residue  is  2TTi ,  we  have  finally 

lim  Ij^  =  2TTi. 

Clearly,   the  integral   I_  convprges  'boimdedly  to  its   limit:   therefore  from  Theorem  VI 
we  know  that  the  limit   in  the  sense  of  distrihutions   of   1^  is   2TTiH(-x). 

V/e  wish  to   investigate  what   is   the  result  of  differentiating   the   limit 
of   Ij^  with  respect  to  x.     From  Theorem  VII  we  have 
-2ni   6(x)=|^     limlj^=limfj     I^  . 

and  from  Theorem  X  we  have 

d        T  r  .    ik3c,, 

-     lH=/s^e       dk; 

consequently,  we  get  the  following  useful  result: 

(2.?8)    lim/g  e^^dk  =-2tt6(x). 

Differentiating  (2.38)  repeatedly  with  respect  to  x,   we  get 

(2.39)    lim/s  (ik)%^''°'dk=-2TT6(^"'^^(x). 

Fere  e^'^^-^-^Jr)  means  the  (n+l)-th  derivative  of  6(x). 

By  Cauchy's  Theorem,  the  integral  over  S  is  the  same  as  the  integral  from 
-R  to  R  along  the  r^al  axis;  therefore  (2.39)  hecomes 

(2.^40)        lim   /  (ik)V^dk=  2tt6(°-'^'^(x). 

R->oo  "^  -R 

Prohlems 

~"  ikx 

2.17  Prove  that  the  limit  in  the  sense  of  distri'butions  of  e    as  k  approaches 

infinity  is  zero.   Prove  also  thpt 

, .   ,  n  ikx   « 
lim  k  e    =0 

k— >oo 
for  n  =  1,2,3,...   Hint.   Consider  y(j(x)e  ^dx  end  use  integration  hy  parts. 

2.18  Prove  the  statement  in  the  text  that  for  0  <  Q  <  "n/z^ 

sin  e  >  -  e. 

—   TT 

2 

Flnt.  Consider  the  gra-^h  of  y  =  sinx  and  y  =  -  i. 
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2.19  Prove  that 


/e^^^dJc  =   lim    j1-^K^^6k  =  it6(x)  +  f 


lim 

E->oo  '  o  €— >0  '  o 

2,20  Use  Problem  2.19  to  evaluate 


/(Hc)'^e 


Il->00       0 

Values  of  DlstrlMtlona  and  Chanre  of  Variables 

It  has  been  stated  previously  that  in  contrast  to  the  case  of  functions  a 
distribution  does  not  have  values  for  different  values  of  x.  Fowever,  for  some  pur- 
poses it  is  convenient  to  define  the  values  of  a  distribution.  First,  we  define  the 
support  of  a  testing  function  0(x)  as  the  closure  of  the  set  of  TPoints  over  which 
t>(j-)  4   0.  We  now  dpfine  the  sense  in  which  a  distribution  ha,g  the  value  zero. 

A  distribution  T  is  zero  In  an  open  set*  -^of  T(^)  =  0  for  every  testing 
function  0  whose  support  is  contained  in  i\-.  Two  distributions  T.  and  !„  are  said 
to  be  equal   in  iLlf  T  -T^  Is  zero  in  JL.   Notice  that  these  definitions  do  not  give 
the  vslup  of  the  distribution  at  8  point  but  only  in  the  nei|P-hborhood  of  a  point.  Also, 
notice  that  these  values  ar?  found  by  using  testing  functions  wnich  are  zero  except 
for  the  neii^hborhood  in  which  we  are  interested. 

We  define  the  support  of  a  distribution  as  the  closed  set  which  is  the 
complement  of  the  largest  open  set  in  which  T  =  0.   If  F  is  the  support  of  T  then  It  If 
clear  that  the  value  of  T(j!5)  will  depend  only  on  the  values  that  ^   takes  in  F.  This 
permits  us  to  extend  the  class  of  testing  fijnctions  v/hen  T  has  a  bounded  support, 
for  then,  since  T(^)  will  not  depend  on  the  values  of  0  for  larre  x,  we  need 
no  longer  require  that  a  testing  function -vanish  outside  a  finite  interval. 
In  the  lang-uage  of  symbolic  functions  we  may  express  this  fsct  by  sayin^  that 
if  T  has  a  bounded  support,  then  t(x)  vanishes  outside  a  finite  xntorvGl  and  so 

T(0)  =/t(x)«5(x)di 
is  really  an  integral  over  a  finite  interval. 

Let  us  investigate  the  valufis  of  6(x).   In  the  neighborhood  of  any  point 
not  the  origin  its  value  is  zero.   In  the  neighborhood  of  the  origin  the  value  of 
6(y)  is  not  defined;  consequently,  the  support  of  P(x)  is  the  origin.   Since  the  sup- 
port is  a  bounded  set,  the  6-functional  may  be  applied  to  any  continuous  function  with- 
out any  restriction  as  to  behaviour  for  large  values  of  x.   Similar  conclusions  apply 
to  any  derivative  of  the  6-functional. 


*  An  open  set  -&~.is  a  set  of  points  such  that  if  P  belongs  to  -A^  then  all  the  points 
of  some  interval  containing  P  belong  to  -S^. 
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We  now  wish  to  define  the  meaning  of  a  change  of  variable  in  a  distribution 
or  in  a  symbolic  function.  Let  f (x)  be  an  absolutely  continuous  monotonlc  function  of 
X  and  such  that  f '(t)  /  0;  then  if  t(x)  is  eji  integrable  function  and  il5(x)  a  testing 
function  we  have 


(2.i4l)      /t(fU))0(x)d^  =/t(y)|«(x)  -p^ 


by  making  the  change  of  vrriable  y  =  f (x).  The  absolute  value- of  the  derivative  is 
needed  in  order  to  ensure  that  the  range  of  integration  should  always  be  from  -  oo  to 

+  00. 

We  use  (2.^1)  to  define  t(f(j-))  in  case  t(z)  is  not  an  integrable  function. 
-) 
For  expmple,  6(x  +?3c-^)  is  defined  by 


/  6(x'\3x-/*)0U)d^  =/6(y)  -^  dy  =  ^ 


3(x2+l) 


since  x  =  1  when  y  =  C.     We  then  have 

In  a  simil-^r  way  we  can  show  thnt 

6(x-xJ 
(2./^2)  6(f(x))  = 


If'^o^l 


where  x  is  the  root,  if  any  exists,  of  the  eaijption 

f(x)  =  0. 
We  shall  extend  (2.^2)  to  the  case  where  f(x)  i'  not  monotonic,  i.e.,  a 
case  when  f '(x)  may  equal  zero.  We  put 

The  numbers  x,  are  the  zeros  of  the  equation 

f (x)  =  0. 
As  an  illustration  of  (2.ii2)  we  quote  a  formula  given  by  Dirac; 

6(x^-a^)  =  ^   [6(x-a)+6(x+a2. 

Finally,  it  is  easy  to  see  from  (2.^1)  that 

|^t(f(x))  =  t'(f(x))f'(x)  , 

which  is  the  same  as  the  usxial  formula  for  differentiating  a  function  of  a  functions. 
It  is  to  ■^e  noted  that  here  the  derivative  is  in  the  sense  of  distributions. 
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Problems 

2.21  Suppose  f (i)  is  an  absolutely  continuous  function;  then  H(f (x))  equals 
1  when  f (x)  >  0  and  eq\ials  zero  when  f (x)  <  0.  Use  the  definition  of  the  de- 
rivative of  a  distribution  to  evaluate  r—  H(f (x))  and  show  that 

ox 

^H(f(x))  =  f(x)6(f(x)), 

where  6(f(x))  is  defined  by  (2.i^2). 

2.22.  Start  with  the  power  series  for  loe(l-Z)  and  show  that 

eo  ^n+l 


f  ^(^1)  =  ^^-^)  ^°^^^-^^  *  "" 


Prove  that  the  series  converges  uniformly  and  boiandedly  for  all  Z  on  the  unit 

cirlce. 

2.23  Prove  that  „ 

00   in9     T                 tt-(«-2tt[-^]) 
n  2-^  =  -  i  loe  (2-2cos«)  +  i  3  ^    . 


Here  L"^jnis^ii8  *^®  largest  integer  in  0/2tt.   (Hint.  Put  Z  =  e   in  Problem 


^Dmeans    the  la ^ '~        '— ^     --   "  -   -^« 

2.22  and  differentiate.) 

2,2^    Prove  that   in  the  sense  of  distributions, 

f  cflsin     =.|i,g(2.2co8©) 

f  smen     .H      ^'^^^ 
r-      n  2  ~         2 

sin  n9  =  -  —  cot  — 
1  /i  ^ 

^  +  ^  cos  ne  =  TT  2_  6(e-2n1c). 

^1  -00 

(Hint.   Use  Problem  2.23  end  Theorem  IX.) 


\' 
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Integration 

Suppose  a  distri'bution  S  is  given  ?jid  we  wish  to  find  a  distri^butioa  T  such 
that 
(Z.l^h)  T'  =  S. 

If  this  were  a  relation  between  functions,  there  would  be  an  infinite  number  of  solutions 
of  (2.^^)  differing  from  each  other  by  a  constant.  We  shall  show  that  the  same  is  trae 
for  distributions  or,  what  is  equivalent,  that  the  only  solution  of 
(2.^5)  T'  =  0 

is  T  equal  to  a  constant. 

The  proof  follows  very  closely  the  proof  in  Chapter  I  that  a  bounded  linear 
functional  is  a  scalar  product.  Suppose  T  satisfies  (2.^-5);  then  consider  the  null 
space  of  T,  i.e.,  the  set  of  all  testing  functions  >\i{x)   such  that  T(\|;)  =  0.  From 

(2M) 

T«(J5)  =  -Tt^')  =  0 
for  any  testing  function  0(x).   This  shows  that  the  set  of  derivatives  of  testing 
functions  belong  to  the  null  space  of  T.   Now,  when  is  \f   equal  to  some  ^'7   If  it  is 
equal,  then 

/  y/  dx  =  /  0'djr  =  0, 
since  i(.x)   vanishes  outside  a  finite  interval. 

Suppose,  now,  th?t  0  (x)  is  a  testing  function  such  that 

/0^(x)dx  =  1. 
Put 

jii{x)d.x  =  a; 
then  t{x)-o0    (x)  is  a  testing  function  such  that 

yO5(x)-a0^(x[ldx  =  0; 

consequently  0(i)-a^  (x)  =  \ffix)   Is  the  derivative  of  some  testing  function  and  there- 
fire  belongs  to  the  null  space  of  T.  ¥e  now  have 

0  =  T(v|/)  =  T()6-a^^)  =  T(0)  -  odr(0^) 
or 


T(0)  =  aT(0^)  =  T(0^)y^dx  . 


This  proves  that  T  equals  the  constant  T(0  ). 
Convolutions  or  Folding  of  Distributions. 

We  have  seen  that  in  general  the  product  of  two  distributions  does  not  exist. 

However,  a  new  type  of  product,  called  the  convolution  or  the  Faltung  or  the  fold*, 
*.  This  term  is  due  to  Professor  Ewald, 
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can  "be   defined.  For  integra'ble  function  f(x)  and  e(x),  the  fold  h(i)  is  given  by  the 
formula 

{2.h6)     h(x)  =/f(i-y)g(y)dy  =/f(y)g(x-y)dy. 
We  write 

h  =  f*g. 

It  is  easy  to  show  that  this  product  is  hoth  commutative  and  associative. 

To  extend  (2.'j6)  to  distributions,  a  functional  definition  must  be  found. 
Let  ^{x)   be  any  testing  function;  then 

h(|«)  =/i«(sc)dx/f(x-y)g(y)dy  =  ^0(xjf  (i-y)g(y)dydi 

where  we  have  put  ^=  x-y,  \=  J ,     The  last  integral  may  be  considered  as  a  repeated 
functional.   Start  with  ^(^+l).  Take  its  g-functional  with  respect  to  \;    the  result 
will  be  a  function  off,  say  ^i{%) .     This  function  is  infinitely  often  different  iable. 
If  g(x)  has  F)  bounded  support,  then  \v{.%)   will  vanish  outside  a  finite  interval  and  so 
it  will  be  a  testing  function.   Now,  by  (2.^7) 

h(0)  =  f  (xp). 
If  g(x)  does  not  have  a  bounded  support,  then  v(5)  need  not  vanish  outside  a  finite 
interval,  so  it  will  not  be  a  testing  function.  However,  f (n;)  will  still  exist  if  f (x) 
has  a  bounded  support.  We  conclude  then  that  h(0)  always  exists  if  either  f (x)  or  g(i) 
has  a  bounded  support. 

Let  S  and  T  be  two  arbitrary  distributions,  one  of  which  has  a  bounded  support, 
The  fold  of  S  and  T  is  defined  as  follows: 
iZ^m)  S*T(0)  =  S^i^T>^[]5(^  +  \a  )  . 

It  is  easy  to  show  that  the  fold  is  both  commutative  and  associative. 
As  an  illustration  of  folding,  consider  6*T.   By  (2.W) 

6*T(0)  =  6^[T^[j!i(?  +  »()]^  =  \[6jii5(f+>^)]]=  T^[^(^)]  ^  T(0)  ; 
therefore 

iZM)  6*T  =  T. 

with  symbolic  functions  this  wcjld  be  written 

t(x)  =/6(i-y)t(y)dy. 
Similarly,  the  fold  6'*T  can  be  evaluated.  By  (2.^8) 

6'*T(0)  =  T^^e^  Q5(<+lI]}  =  -T^B'Wa  =  T-(0)  ; 
therefore 
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6'*T  =  T>. 
A  aimilpr  ar^-ument  shows  tJnet  if  D  is  any  differential  operator,  then 

(2.50)  D6*T  =  DT. 

Also,  by  (2.50)  and  the  associative  law, 

(2.51)  D(S*T)  =  D6*(S*T)  =  (D5*S)*T  =  DS«T  =  S»DT; 

thus  we  see  that  to  differentiate  a  fold,  only  one  of  the  factors  need  be  differentiated. 

These  relations  are  very  useful  in  the  theory  of  differential  equations.  The 
solution  of  the  differential  equation 

(2.52)  DT  =  E, 

where  R  is  a  given  distribution,  can  be  found  once  the  solution  of  the  following 
equation: 

(2.53)  DS  =  6 
is  known.  For  if  we  put 

T  =  S*R 
then,  by  (2.51),  (2.52)^ and  (2.^9)  we  find  that 

"DT  =  DS*R  =  6*R  =  R. 
The  distribution  S  satisfying  (2.53)  will  be  called  the  fundamental  solution 
of  the  differential  equation  (2.52). 


Problems 

ow  that  for  n  >  2,  S  =  - 

(n-2)x' 


2.25  Shov;  that  for  n  >  2,  S  =  -  ^   '  is  a  fundamental  solution  of  the  equation 


(I'^-^f)'  =  r. 

"For  n  =  2,  S  =  E(x)lo/?  x. 

Di3tribution  in  More  than  One  Dimension 

All  the  previous  ideas  can  be  carried  over  to  an  n-dimensional  space.  A 
testing  function  ^  is  defined  to  be  a  function  of  the  n  variables  x^ ,  x-,  ...,i  which 
is  Infinitely  different lable  and  which  vanishes  outside  some  bounded  region.  A  dis- 
tribution is  a  linear  continuous  functional  on  the  spece  of  testinc  functions.  Any 
function  f  (x-.  ,x  , . . .  ,x  )  which  is  stmimrble  in  n  dimensions  defines  a  distribution  by 
the  following  fonnula: 


f(^)  =f,..ff{x^ x^)  0(Xj^ x^)  dx^ dx^. 


The  derivative  of  a  distribution  T  is  defined  by  the  following  formula  which  is  slmi- 
Ipr  to  that  used  in  the  case  of  distributions  of  one  variable: 

All  the  theorems  previously  stated  for  distribution  in  one  dimensions  also  hold  for  dis- 
tributions in  n  dimensions. 
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Besides  the  distributions  we  have  introduced  above,  we  shall  also  consider 
the  direct  product  of  a  distribution  in  one  space  and  a  distribution  in  another  space. 
In  order  to  explain  this  concept,  we  shall  work  with  a  simple  illustration. 

Suppose  th,^  t  S  is  a  distribution  on  the  one-dimensional  space  of  testing 
functions  0(x)  while  T  is  a  distribution  on  the  one-dimensional  space  of  testing  func- 
tions  w(y).  From  the  one-dimensional  testing  functions  ^{x)   and  'i;(y)  we  may  form  two- 
dimensional  testing  functions  9(x,y)  by  the  following  device: 

i=l 

It   is  easy  to   shoi^  that  every  two-dimensional   testing  function  of  x  and  y  can  be  approx- 
imated by  sums  such  as   thpt   in    (2.5^0. 

We  now  show  thi?t  we  may  define  a   two-dimensional  distribution  W       on  the  space 

xy 

of  two-dimensional  testing  functions  6(x,y).  This  distribution  has  the  property  that 

^^'""^^  V  05(3tV(y[!  =  s^(0)Ty(v!/). 

Consider  any  two-dimensional  testing  fijnction  ©(x,y)  as  a.  function  of  x  depending  on  a 
parameter  y.  For  fixed  y,  Q(x,y)  is  a  testing  function  of  x  and  then 

is  a  function  of  y.  From  Theorem  X  it  follows  that  *(y)  is  a  testing  function  in  y  and 
we  may  then  define 


Ty[J(ya  =  Ty[s^[e(x.ya]  =  w^^y(e) 


It   is  clear  that 

w[ii(x>v(ya  =  Ty[s^I]5(x)M/(yi]]=  Ty&s^(^|]  =  Ty(w)s^()«). 

The  above  construction  shows  thpt  there  exists  a  distribution  W  having  the 
property  (2.55).  V^e  call  W  tne  direct  product  of  S  and  T.  This  distriViti^r.  must  '^e 
unique  since  it  has  a  unique  value  on  every  sum  of  the  form  (2.5^). 

Note  also  that  W  could  just  as  well  be  defined  by  the  formula: 
W[e(x,y[!  =  S^[Ty(^(x.y[l. 

Since  W  is  unique,  this  imples  that 

♦  1,  Schwartz,  Theorle  des  Platributions.  op.cit. ,  tome  1,  p. 107. 
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s,[TyCe(^.y)3]  =  Ty[s^[e(z,ya] 


or,   expressed  in  terms  of  57111110110  functions, 

(2.56)  /s(x)dx/t(y)e(-,y)dj-  = /t  (y)dy/s(x)e(3r,y)<ix. 

We  may  s^ummarize  o^xr   results  in 
Theorem  XI.   If  S  is  a  distril^-ation  over  the  space  of  testing  functions  0(t:)  and  T 
Is  a  distri'bution  over  the  sjiace  of  testing  functions  H/(y),  then  there  exists  a  unique 
di3tri"hution  '■''  over  the  space  of  two-dimensional  testing  functions  Q(x.y)  with  the 
property  that 

w|0(^V(y)l  =  S(?5)  T(,j,)  . 

wCeCx.yfl  =  S^[Ty[i.(x.ya]  =  Ty[s^[:o(x.ylf|  . 

Note  that  the  last  result  or  what  is  eciuivalent,  formula  (2. ^"7),  amounts  to 
a  statement  of  the  fact  that  we  may  interchange  the  order  of  integration  of  symholic 
functions. 

The  concept  of  direct  product  of  dlstrlhutlons  can  ohviouely  be  extended  to 
the  case  where  we  start  with  distributions  defined  on  m-  and  n-dimensional  space,  and 
to  the  case  where  we  have  more  than  two  factors.  We  leave  the  details  of  these  cases 
to  the  reader. 
S-Timctions  in  Different  Coordinate  Systems 

The  n-dimensional  6-function  &{x,,x^,...,x^)   is  the  symbolic  function  for  the 
distribution  which  assigns  the  value  i!5(0,0, . . .  ,0)  to  the  testing  function  (^(x^.Xg, . . .  ,x^). 
Since  the  n-dimensional  6-function  applied  to  the  product  of  one  dimensional  testing 
functions  i^{T^)^^{T^)    '-'^^(\)   gives  the  value  f)^{Q)p^iO). .  .f)^{0) ,   we  see  from  Theorem 
XI  that 

(2.57)  eCx^^.Xg x^)  =  6(x^)6(x2)...6(x^). 

Consider  now  the  special  case  of  two  dimensional  space  and  suppose  that  we 
change  coordinates  from  ^..x  to  ^^,  S  ^  usinF  the  formulas: 

Here  we  assume  that  u  and  v  are  continuous  differentiable  one-valued  functions  in  the 
ordinary  sense  and  we  also  assume  that  ^i  =  Pn»  ^o  ~  ^2  '^°^^®^P°'^'^^  *°  ^1~^*^2~'^2* 
The  integral 

by  the  change  of  coordinates,  becomes 
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(,.,8)  //e[«(f,.|p-={]«['(fr?2>-<'2]*''-'"'''^^'^^  ''''^••^*' 

where  J,  the  Jacobian  of  the  coordinate  transformation,  is  ^iven  by  the  formula 
Equation  (2.58)  "hows  that  the  symbolic  function 

assigns  to  any  testing  function  ^ -,*  5p  »  *^®  value  of  that  function  at  the  points 
where  u  =  cl.  ,  v=  a^,  or,  what  is  the  same,  where  5-  "  ^i  »  5"?  ~  P?»  consequently,  we 
may  write 

or,  finally,  if  J  ^^  0 

6(VPj6(f2-p2) 
(2. 59)     6{x^-cx^)t(x^-rt^)  =   \j\  • 

As  an  illustration  of  this  theorem  consider  the  transformation  from  rectangular 
coordinates  x,y  to  polrr  coordinates  r,e.   If  0(x,y)  is  a  testing  function,  then 


// 


6(x^x^)6(y-y  )0(x,y)dxdy  =  0(x  ,y  ). 

0        0  0   0 


Suppose   that   the  point  with  rectan^mlsr  coordins-tes  x   ,   y     hss  polar  coordinates   r   ,   Q 

o        0  0        0 

and  sup-nose  that  ?(r,0)   is   0(x,y)  expressed  in  polar  coordinates.     Since 


// 


6(r-r   )6(Q-e   ) 

2_ 2_    ^(r,e)rdrde  =  ^(r^.e^). 


and  since  ?5(r    ,Q   )  =  0(z    ,y   ),    it   is  clear  that 
0       o  0       0 

6(r-r   )6(e-9   ) 
(2.60)  6(x-x^)6(y-y^)  =  2_ 2_      . 

This  is  exactly  the  result  that  would  be  obtained  hy  an  application  of  formula  (2.59). 

We  now  investigate  what  happens  to  formula  (2.59)  when  J  =  0.   In  such  a  case 
the  transformation  from  X-.Xp  to  Tit  ^p  breaks  down  and  is  no  longer  one-to-one.  For 
example,  in  the  transformation  from  rectangular  coordinates  to  polar  coordinates,  the 
Jacobian  becomes  zero  at  the  origin  x  =  y  =  0  and  we  know  that  there  the  transformation 
is  no  longer  one-valued  since  r  =  0  and  Q  may  have  any  value.  This  situation  is  typical 
of  the  general  case.  We  shall  call  a  coordinate  such  as  Q,   which  is  either  aany-val-.ied 
of  else  has  no  determinate  value  at  a  singular  point  of  the  transformation,  an  ignorable 
coordinate. 

Suppose  that  x,=a- ,  Xp=ap  is  a  singular  point  of  the  transformation^ that ^ „  1b 
an  ignorable  coordinate  there  and  suppose  that  the  testing  function  0(x, ,x  )  in.  x^,X2 
space  becomes  the  testing  function  ^(^.fg)  ^^   fi  fp   ^P^°®'  ^^^^   ^1~°^*  ^2~°'2  ^°^^^^~ 
ponds  to  f  ^  =  &^,  ^2  indeterminate,  and  also  f)ia,^,a^)   is  a  function  ^(P^)  only  of  ?^. 
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The  sii/mbolic  function  corresponding  to  6(3:^-0- ) 6 (x-ttg)  wili  con3e(iuently  "be   a  symbolic 
function  t(^.  )  only  of  ^,.  Formula  (2.58)  then  becomes 

(2.61)    0(0^,02)  =yy  6(u-a^)6(v-a2)0(u,v)  \j\    d  ?  ^cl  ^  ^ 

Put 

/  lj|dJ2  =  J^; 

then  (2.60)  will  he  satisfied  if  we  put 

6(5  tP,) 


consequently,  in  esse  J  =  0,  when  x,   =   a-,  x^  =   a-,  we  have 

6(V^ 
(2.62)         tix^-a^)6{x^-ci^)   =  j- 


6(VP^) 


'1 

For  example,  in  the  transformation  from  rectangular  to  polar  coordinstes,  Q   is  an  ig- 

norahle  coordinate  at  the  origin,  and 
2tt 
/    rdQ  =  2nr; 

consequently, 

6(x)6(y)  =  ^   . 

Similar  results  hold  for  transformations  in  n-dimensional  space.  We  state 

then  in 

Theorem  XII.  Let  t, ,x   , . . . ,x     "be  n-dimensiooal  rectangular  coordinates  and  let 

€it  5'o...«.  ?  be  any  other  coordinate  system,  not  necess.irily  rectangular,  with 
-i_I C-^  r  n  ^ — ■ — ■ — 

n-dimenslonal  volume  element  JdS^d?  .  ..dg  .   Suppose  that  the  point  F  with  coordinates 

r,  =  rr     X  =a          .  ,x  =a.     has  coordinates  f  ■,  =  P., ,. ..  ,$  =P  and  that  J  -^  0  at  P;  then 
— 1 —  1—  2 — r n — ^n  1 1-* ^  n — n  ' 

6(r-p  )...5(^-p  ) 

(2.63)  B(x^-a^)    ...  Hx^-a.J   =  ^--^--j ^^-^  . 

If  J  =  0  8t  the  point  P,  suppose  that  the  eqtiations  >   =  3^,...,^  =P^  define  P  and 
therefore  the  coordinates  ^v^.Tf-Sj,  ^^e  ignorable.   Put 

so  that  J  is  the  integral  of  the  Jacobian  over  the  Ignorable  coordinates;  then 


6(?^-p^)...6(rj^-ej^) 


(2.6')     6(x^-a^)...6(x^-a,^)  =  , 

k 

The  proof  of  this  theorem  will  be  left  to  the  reader.  We  shall  illustrate 

it  by  considering  the  transformation  from  three-dimensional  rectangular  eoordinptes 

x,y,z  to  spherical  coordinates  r,9,j!5.   Here  J  =  r^sinO  which  vanishes  for  points  on 
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on  the  poler  axis  where  0  is  ignora"hle  and  also  at  the  origin  where  both  Q   and  0  are 

t 

0 


ipnorable.      If  p  has  coordinates  x   ,    y   ,    z     in  rectan£"alar  coordinr.tes  and   t   ,   Q    .    i 

0         0         0  0         o         c 

In  spherical  coordinates,   then 


6(r-r   )6(e-e   )6(0-^   ) 

6(x-xj6(y-y   )6(.-z    )  =  ^ ^ ^    .  when  r^j^O.  ^0 

°  °  °  r'^  sine 

6(r-r^)6(e) 

(2.65)  =  Zm •  ^^^'^    W°'   \^^  °'*  'o''°' V°' 

Zttt  sino 


Problems 


2      »  wh«n  I  =y  =z  =0  or  G    =  r     =0, 

UttT  0*00  0  o 


2.26  Verify  the  relations   in  eq-uation    (2.65)  hy  operating  with  both  sides   in  testing 
functions. 

2.27  Show  that  in  two  dimensions 

6  (r-r  )  .  2tt 

— ^-    =  /         6(x-r  cos©)6(y-r  sine)de   . 

r  »/    Q  0  o 

Hint.  Integrate  equation  (2.60)  with  respect  to  9  . 

2.28  Verify  formula  (2.65)  by  ap  lying  it  to  a  testing  function  which  is  equal  to 

-x-y-z  .   4..    .     2^  2^  2   , 

e      in  the  sphere  x  +y  +z  =  1. 

2.29  Verify  Theorem  XII  by  applying  each  side  of  (2.63)  and  (2.6^)  to  a  testing 
function. 

Courier  Integrals 

The  Fourier  integral  for  a  distribution  T  is  defined  to  be  T  (e    )  if  it 

-ikx 
exists.  The  theory  is  complicated  becatise  of  the  fact  that  e     is  not  a  testing 

—ikx 
function.   However,  we  can  approximate  e     by  a  sequence  of  testin/^  functions  0  (k,x) 

which  vanish  outside  an  interval  whose  len<p-th  goes  to  infinity  with  n,  and  then  we  define 

{?.66)  T^(e-^^)  =  lim  T^(0^). 

where  the  limit  is  in  the  sense  of  distributions.   Of  course  if  T  has  a  compact  support, 

the  limit  in  (2.66)  and  therefore  also  the  Pourier  integral  always  exist.  For  example, 

if  T  is  the  6-function  its  Fourier  integral  is 

/6(k)e"^^dk  =  1. 
If  1  is  considered  as  defining  a  distribution,  its  Fourier  integral  will 
be  defined  by  (2.66).  The  right  side  of  (2.5^+)  is  essentially 


This  limit  may  be  obtained  by  using  the  theory  of  the  Dirichlet  kernel*  but  we  see  that 


lim       /       e^^dk. 

^     i/  -n 
n— >Q0 


•  See  Ti'tchmarsh,   Theory  of  Functions    op*cit.  , 
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it  is  merely  the  case  n  =  0  of  formula  (2.39).  We  have  then 

(2.67)  fe'^^dk  =   2tt6(x) 
or,  replacing  x  "by  x-x',  we  get 

y.^-ik(x-x')^^  2tt&(x-x'). 

Consider  an  arbitrary  di8triT)ution  T  and  take  the  fold  of  (2.6?)  with  T.  Using 
symholic  functions,  we  have 

/t(x)dx/e-^^^''-^' V  =  2i/t(x)6(x-x')dx  =  2TTt(x«); 

hut  with  the  help  of  Theorem  XI  we  may  interchange  the  order  of  integration  and 
find  that 

(2.68)  2TTt(x')  =/e^^'dky  t(x)e~^'^dx. 

The  inner  integral  is  exactly  the  Fourier  integral  of  the  distribution  T^  as  we 
have  defined  it  in  (2.6  6),  hut  now  we  consider  T  as  acting  on  functions  of  i  and 
not  on  functions  of  k. 
Put 

(2.69)  T(k)  =/  t(x)e"^^dx 

where,  as  always,  the  right  side  is  to  be  understood  in  the  sense  indicated  in 
(2.66).  We  shall  call  T(k)  the  Fourier  transform  of  t(x).  Formula  (2,^)  now 
gives  an  inversion  theorem  for  the  Fourier  transform,  namely 

(2.70)  t(x)  =  ^  /T(k)e^'^dk. 

Fote  this  inversion  formula  is  the  Fourier  transform  of  T(k)  except  for  the  change 
of  sign  of  the  imaginary  unit  and  the  factor  1/2tt. 

We  shall  illustrate  the  inversion  formula  by  considering-  the  Fourier 
integral  given  in  (2.39).   It  states  that 

X  e    dx  =  2TTi  6   (k)^ 

consequently  the  Fourier  transform  of  x  is  2ni  6   (k). 

Now,  the  inversion  formula  would  give  the  following: 

n   1  /*  -._jn.(n)/,  \   ikx,, 
X  =  ^zj   2TTi  6^  '(k)  e   dk. 

This  formula  is  clesrly  correct  since 

/e^'^^k)  J5(k)dk  =  (-)'^^^^^o). 

It  should  also  be  noticed  that  as  a  consequence  of  (2.39)  every  polynomial 
has  a  Fourier  transform  and  that  this  transform  is  a  linear  combination  of  deriva- 
tives of  the  delta  functions.  We  state  without  proof  th^-^it  the  fact  that  every  inte- 
grable  function  f (x)  which  satisfies  the  inequality 
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for  all  c  >  0  and  all  x   sufficiently  large  in  absolute  value,  has  a  Fourier  trans- 
form. This  Fourier  transform  will  in  general  be  a  distribution.  We  shall  not  dis- 
cuss under  what  conditions  a  distribution  has  a  Fourier  transform.   In  applications 
it  is  usually  not  difficult  to  decide  whether  the  Fourier  transform  exists.   If  it 
does  exist,  the  inversion  formulas  (2.69)  end  (2.70)  will  hold. 
Problem: 


2.30  Prove  that        -       .,        ^. 
■^  f  -itac.      2i 

J   sgnx  e    dx  =  -  7j^  J 

and  prove  also  the  inverse  formula 


2tt/ 


-2ie   dk 

=  sgn  X. 


2tt/      <k 
Properties  of  the  Fourier  Transform 

We  assijme  that  all  distributions  considered  In  this  section  have  Fourier 
transforms.  The  first  property  we  sha]l  discuss  is  fiven  in 

TheoT-em  XIII.  If  T  is  a  distribution  with  compact  support,  then  the  Fourier  trans- 
form of  T  is  an  entire  analytic  function. 

Suppose  that  the  support  of  T  is  contained  in  the  interval  -R  <  x  <  H.  Let 
0.  (x)  be  any  testing  function  which  is  identically  equal  to  e~  '   for  -R  <  x  <  E;  then 
the  Fourier  transform  of  T  is  given  by  the  following  formulp  which  is  valid  for  all 
real  and  complex  k. 

rCk)  =  T(e-^^)  =  T(0^). 
By  Theorem  s^   the  derivative  of  ^(k)  is 

r(k)  =  f^T(^,p  =  T(^0^)  =  -lT(xe-^^) 

which  again  exists  for  all  real  and  complex  k.  Since  the  derivative  of  '^(k)  exists 
and  is  unique  for  all  conrplex  vslu»=s  of  k,  it  is  proven  that  M^k)  is  an  entire  analy- 
tic function  of  k. 

Note  that  since 

T(x0)  =  xT((i), 
we  have  also  proved  that  the  '''ourler  transform  of  xT  is  1  'if'(k).  From  this  it  follows 
easily  that  the  Fourier  transform  of  x  T  is  1  X        (k).   Conversely,  it  can  be  shown 
ti^at  the  Fourier  trsneform  of  T    is  inc  "^(k).  We  re-stste  these  properties  of  the 
■^ourler  transform  as 

Theorem  XIV.  If  T  is  a  distribution  which  has  a  Fourier  transform  t^(k).  then  x°T  has 
Fourier  transform  l'^^^"^(k)  and  T^'^^  has  Fourier  transform  i^k*^  t(k). 

The  proof  above  used  the  fact  that  T  has  a  compact  support,  but  it  is  easy 
to  show  that  this  restriction  is  unnecessary. 
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The  last  property  of  the  To-urier  integral  which  we  shall  discuss  is  the 
convolution  or  folding  property.  We  have 

Theorem  XV.   Suppose  thpt  T  and  S  are  any  two  distrihutions  one  of  which  has  a  com- 
pset  support,  and- suppose  that  the  Fourier  trgnsform  of  T  is  ir(k)  and  the  Fourier 
transform  of  S  is  <r(k);  then  the  Fourier  transform  of  the  fold  of  T  and  S  is  the 

product  ^(k)q-(k). 

?Tote  that  the  product  T(k)  Q-(k)  always  makes  sense  since  hy  Theorem  XIII 
and  the  hypothesis  one  of  the  functions  is  analytic. 

The  proof  of  Theorem  XV  is  easy.   By  definition  the  fold  of  T  and  S  is 
T*S((!»  =  T^[Sy05(x+yi]:] 

and  hy  definition  the  Fourier  transform  of-  the  fold  is 

T*s(e-^^)  =  T  CS  (e-^'^-^'^a  =  T  (e-^^)S^(e-^^)  =  r(k)  <r(k). 

X    y  A  Jr 

This  proves  the  theorem. 
Pro'blems : 

2.J1     Find  the  Fourier  transform  of  \y.\  .     Use  thp  result  of  Problem  2.30  to  show  that 

the  result  is  i  times  the  derivative  of  the  Fourier  transform  of  sgn  x.   Compare  this 

with  Theorem  XIV. 

2.32  Prove  Theorem  XIV  for  an  arbitrary  distribution  that  has  a  Fourier  transform. 

P.."^"}     Prove  the  converse  of  Theorem  XV:  If  the  product  of  the  distributions  T  and  S 

has  a  meaning  and  if  T  and  S  have  Fourier  transforms  t  (k)  and  C (k)  respectively, 

then  the  Fourier  transform  of  the  product  TS  is  the  fold  of  t (k)  and  (r(k). 

?,3h     Use  the  result  of  Problem  2.?^   and  formula  (2.50)  to  prove  Theorem  XIV  again. 

Fpijirier  Transforms  in  Three  Dimensions 

The  Fourier  integial  and  the  inversion  theorem  can  be  extended  to  distri- 
butions in  n-dl!!iensions.  V/e  use  the  n-dimensl6nal  equivalent  of  (2.47).  namely 

/...yexpG-i(k^x^+...k^i^3<ik;iL*"'^n  ""  (2TT)*^f-(_x^). .  .6(x^). 
Replace  x.  by  x.-^.  for  1  =  1,2,..  ,n  and  wS  get 

y...yexp[-i^kj^(x^-f^)+...+kj^(x^-£^)J]dk^...dkj^  =  (2n)^6(ij^-£^)...6(x^-^). 

This  formula  can  be  used  in  the  same  way  as  the  corresponding  formula  m  one  dimension 
to  prove  that  if  t(x, ,X-, . . . ,x  )  is  an  n-dimensional  function  (symbolic  or  not)  and  if 
we  define  the  Fourier  transform  by 

(2.71)    1:'(k^,...,k^)=y.  ../t(x^ x^)expG-i(k^x^+...+k^x^2di^...dx^  , 

then  the  following  inversion  theorem  holds: 
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(2. "72)   t(x^,..  .,x^)=(2TT)"y..^t'(ki...  .,kj^)ezp[i(k^x^+...+k^x^[lcLk^...dk^. 

The  three  dimensional  case  is  of  particular  Interest.   Pat  n=3,  let 
i^,x_,x_  he  the  components  of  a  vector  r  and  let  k^,kpjk-  he  the  components  of 
a  vector  k;  tnen  (2.71)  and  (2.72)  hecome 

(2.73)    tr(k)  =y  t(r)e-^^dr 

t(r)  =  (2n)-y  (k)e^^'idk 

Here  dr  and  dk  are  the  volume  elements  in  three-dimensional  apace  and  k.r  indicates 
the  scalar  product  of  the  vectors  k  and  r. 

The  properties  of  the  Fourier  transform  discussed  in  the  previous  section 
can  he  suitshly  generalized  so  that  they  are  valid  for  the  transforms  in  (2.73)-  We 
shall  show  how  this  can  he  done  for  the  convolution  or  folding  property.   Let  t(r) 
and  s(r)  he  f-ijnctions  which  have  Fourier  transforms  "flk)  and  ®~(k)  respectively. 
The  fold  of  t (r)  and  8(r)  is 
yt(r«)9(r-r')dr'. 

The  Fourier  transform  of  the  fold  is 

ye-^^-^r/t(r')3(r-r«)dr-  = /e'^^L-l't  (^ .  )ir/e-^^- (^-^' )s(r-rt  )dr 


=  ye-^=-^'t(r')drye-^^-4(£)d5  =  t<k)cr<k); 


consequently,  we  have  the  ss?ne  result  ss  in  the  one-dimensional  case,  namely,  the 
Fourier  transform  of  a  fold  is  the  product  of  the  Fourier  transforms. 

In  many  applications  t(r)  is  a  function  solely  of  the  length  r  of  the 
vector  r.   In  that  case  the  formula  (2.73)  ^o^  *^6  Fourier  transform  will  he  simpli- 
fied hy  the  introduction  of  spherical  polar  coordinates  r,Q,^.  The  angles  Q,  j) 
are  defined  relative  to  ?  polar  axis  which,  for  simplicity,  we  take  to  be  along  the 

direction  of  the  vector  k.  Formula  (2.73)  now  becomes 

^QD   tt/2  2tt      .,     -  „ 
"   y  '      ^  -ikrcose  2 


t(k)  =/^    /^      J^  t(r)e-^''^'°°"®rSinedrd9d0 


where  k  is   the  length  of  the  vector  k.      If  the   integration  with  respect  to  9  and 
f)  is  carried  out,  we  obtain 


(2.7^)        r(k)  =r(k)  =  hrrf^  t(r)  Slgl  Ar. 
Similarly,  we  obtain  the   inverse  formula: 
(2.75)         b(r)  =  ^   //f(k)  ^    k^dk. 
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Problems 

2,  ?5       T'ind  the  three-dimensional  ?oiirier   transform  of  r     where  r  =    (x  +y  +7    )   '     . 

Hint.    SeTiarate  Problem   (2.20)   into  rpal  and  imaginary  parts    to  obtain  the  required 

integral . 

2.35   Show  that  the  integral  over  six  dimensions 

/  f  f(T^)g(\T^-r^\ )f(r2)dr^dr2 

=  /  !a,(k)!^e(k)dk 

where  n.(k),  P(k)  are  the  Fo^irier  transforms  of  f(r)  and  g(r)  respectively.   Hint. 
Replace  r(r^-rp)  be  an  integral  in  term-,  of  its  Fourier  transform  3(k)  as  follows: 

ed^-r^)  =/e^^->^-^)p(k)dk. 

2.37   Use  Problem  2.36  to  evaluate  the  following  six-dimensional  integral: 

/ye-""!  e-^2  Ir^-r^ldr^dr^  . 

Evaluating  '"'inite  Part  Integrals 

The  concept  of  the  finite  part  of  integrals  which  we  introduced  previously 
(see  p^re  80)  plays  a  very  useful  role  in  the  applications  of  Fourier  integrals.   Be- 
fore discussing  these  applications  we  shall  show  how  in  certain  cases  the  finite  part 
of  integrals  may  be  evaluated  by  means  of  complex  integration. 

Consider  the  finite  part  of  the  integral 
00  -ikx. 


-/ 


dk 


-00   Ckr 

where  n  is  a  positive  integer.   I  is  the  limit  as  E  approaches  infinity  of  the  integral 

B  -ikx 

Let  us   investigate  the  value  of  the   integral  of  over  the   contour    L^  in  the 

complex  k-plene   illustrated  in  Figure  1,     The  contour  C  is  made  up  of   two  semi- 


Figure  1 
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circl'=9,  one  of  radius,  R,  the  other  of  rsdius  c,  and  two  seg-ments  of  the  reel  k- 
axis,  namely  from  -R  to  -e  and  from  e  to  R, 

Since  the  inte,?rand  is  rfpular  inside  C-,  then  by  Cauchy's  theorem  the  integral 
over  C  will  vanjsh.  If  x  is  negative,  the  inteprnl  over  the  semi-circle  of  radius 
R  will  vanish  as  R  approaches  infinity}  we  have  then 


(2.'76) 


-ikx       d     -ikx 
dk  +  / dk. 


The  last  term  is  the  integral  over  the  small  semi-circle  from  -g  to  e. 

Notice  thfit  the  integrals  in  (2.76)  are  ordinary  integrals  and  not  only  the 
finite  parts.   Of  course,  as  e  approach-^s  zero  each  term  in  (2.76)  becomes  in- 
finite^ but  their  sum  remains  7ero.  ''/e  conclude  then  that  the  finite  part  of  the 
first  t^rm  is  equal  to  the  negative  of  the  finite  part  of  the  second;  consequently, 

-ikx. 


/I    -IKX,, 


s  Ck)' 

This  last  inteftral  may  be  evaluated  by  expanding  the  numerator  in  powers  of  x. 

Ve  h.^ve 

^  ^    ^  k^dk 


I  =  _  ^  (.ix)J  / 


Now,  as  e  approaches  zero  consider  the  limit  of  the  integral 


/ 


S  ('k)"" 


dk 


Wien  J  is  less  than  n-1,  the  limit  is  infiiity;  when  j  =  n-1,  the  limit  is  -ni;  when 
j  is  greater  tha.n  n-1,  the  limit  is  jero;  consequently, 

I  =  TTi  (-ir)''"-   • 

for  X  negptive. 

If  X  is  i->ositive,   we  use  a   similar  contour  wit^^   s'emi-clrcles   in  the   lower 
hs-lf-plane  and  we  find   that 

_  .  (    .    vn-1 

I       =    -TTI(-IX) 

Ve  may  write  both  results  in  the  form 

(2.76a)       I  =  iTx''-^(-i)''  sgn  x. 
As  another  ill\istration  consider 
.00  -ilex. 


/ 


dk 
> 


O    I  (T^.V^) 


ic„  >  0 

0 


Replace  the  intep-ral  alonr  the  real  axis  by  the  integral  over  the  contour  C,      illus- 
trated in  T'igure  2. 
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Figure  2. 
Since  there  are  no  Bingularities  Inside  ^,  the  integral  over  Cis  zero. 
For  X  negative,    the  integral  over  the  large  semi-circle  vanishes.  The  integral  over 
the  small  semi-circle  is  -Tii  e  "^^  ;  consequently,  for  x  <  0  we  have 


/ 


00  -ikx,, 
e    dk 

'(k-k  ) 
0   ^   o' 


t+ik 


The   last   term  is   the  contribution  from  the   integral  over  the  imaginpry  axis  where 
we  have  put  k  =  it.     For  x  positive  we  use  a  similar  contour  in  the  lower  half- 
plane  and  we  ohtain 


00     -Ikx 


/    %: 


_dk 


=  -ni  e 


-ik„x 


.00  -tx 


o^       /    £ ° 

~J      t-ik 


dt 


Similar  methods  can  he  applied  to 
.00 


c  oskxdk 
0  '(k^-k^) 


V^e  first  write  the  integrals  in  terms  of  exponentinls  as  follows: 


j  =  i 


=  ^/„- 


00     ikx^  -ikx 
e       +e 


(.^-.^ 


dk 


=  1/2 


00     ikx 
e 


(k^-k^) 


dk 


o   Hk  -k  )  '   0 

This  last  integral  is  replaced  hy  a  contour  integral  indented  at  the  singularities 
and  closed  in  the  upper  or  lower  half -plane  according  as  x  is  positive  or  negative, 
-vsluating  the  contrihution's  at  the  singularities  we  find  that 

J  =  --^  9in(kjx!). 
o 


As  a  final  lllustrntion  consider  the  integral 
K 


'/- 


00  ^-Ikx,, 
e    dk 

00 


Ck)' 


where  u  is  an  ar'Mtrf'ry  positive  member.   In  the  beginning  of  this  section  we  have 
evalimted  the  Inte/Tal  when  i>  is  an  integer.   If  v   is  not  an  integer  then  the  inte- 
gral K  will  not  he  completely  defined  until  we  specify  which  hrancl:  of  the  multi- 
valued function, k  we  are  using  in  the  integral.  For  example,  if  i?  =  1/2  then  when 


k  is  positive  k 


1/? 


may  be  taken  rs  the  positive  square  root;  but  when  k  is  negatlve^does 
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1/2 

k         mean  the  positive  imp^lnary  or  the  negative  ime^insry  sqimre  root?     Sup-^ose 

1/2 

we  specify  thpt  k    should  mean  the  positive  imaginary  sqtuare  root.   Then  we  may 

evaluate  F.,  #„  'by  the  previous  methods  if  we  take  into  account  the  fact  that,  'heceuse 

'  1/2 

of  oior  specification  of  the  value  of  k   ,  the  lower  half  of  the  complex  k-plane 

must  contain  a  branch  cut.  Vfhen  y  is  negative,  the  contour  may  he  closed  in  the 

upper  half-pl8.ne  thu?  avoiding-  the  branch  cut^and  we  find  K^ /p  =  0.  When  x   is 

positive,  the  appropriate  contour  is  illustrated  in  Pifure  3-  The  shaded  line 

along  the  negative  imaginary  axis  represents  the  branch  cut.  The  integrals  over 

the  semi-circles  vanish  but  the  integrals  over  the  branch  cuts  do  not  vanish. 


Figure  3 

The  integral  on  the  left  side  of  the  cut  is  equal,  except  for  a  factor  -1  to  the 

integral  on  the  right  side  of  the  cut.  '•Replace  k  by  -it  in  the  branch  cut  integrsls 

and  we  have  . 

_TTi/ii   /°°elZdt       ^r    -W^ 

0  u 
for  X  >  0. 


1/3 


A  similpr  treatment  must  be  used  for  other  values  of  u.   It  may  seem  'onnecesspry 

i/t 

when  V  =  l/''  say,  for  then  it  woiild  seem  thpt  k    is  positive  when  k  is  positive 
and  negative  when  k  is  negative.   Fowever,  with  this  specification  the  function  k 

is  no  longer  continuous  in  the  com-nlex  k-plane.  This  may  be  seen  by  following  the 

1/? 
values  of  k    along  the  unit  circle  k  =  +1  to  k=  -1.   If  we  start  at  k  =  1  with 

1/3  1/3 

k  '   =  1,  then  if  we  go  along  the  unit  circle  in  the  upper  half-plane  we  get  (-1)  ' 

=6  '-  but  if  we  go  along  the  unit  circle  in  the  lower  half-plane  we  get  (-1)  ' 
=  e~  '  .   In  neither  case  can  we  get  (-1)  '^  =  -1, 


plane 
Since 


1/3 
Nevertheless,  despite  the  fact  that  k  '  '  is  not  continuous  in  the  complex  k- 

,  the  value  of  K^ /_  may  be  obtained  but  we  must  proceed  slightly  differently. 

1/3  '^ 

k    is  an  odd  function  of  k,  we  have 


/ 


e         dk 


-  00 


jrr 


sinkx 


173 


dk 


=/ 


00     -ikx     ikx 
e         -c 


Jh 


dk. 
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1/3 

Let  us  put  the  'branch  cut  of  k    alon^  the  negative  real  ails;  then  one  term  of 

the  inteisrral  may  be  rotated  in  the  u-^per  half-plane  until  it  lies  elonp  the  upper 
ima^iTiBry  axis  while  the  other  is  rotrted  in  the  lower  helf-plane  until  it  also  lies 
plon^  the  lower  ima,?insry  axis.   In  thl»  way  we  find  that  if  x  >  0. 
.-ioo  -ikx  .     loo  itac 

Put  k  =  -it  and  k  =  it  respectively  in  these  integrals  and  we  have 

ttI   ttI       , 
t  '~      Tn  /*  e-^dk 

V3=  ^'     -'    V,   ^iJT   • 

The  last  integral  can  he  found  from  tables  of  the  Ganuna  Function. 

Trom  these  examples  we  may  derive  the  following  method  for  finding  the 
finite  part  of  an  integral: 

Me thod .  Replace  the  integral  by  a  contour  integral  indented  at  the  singularities 
and  closed  in  the  upper  or  lower  half-plane^and  evaluate  the  conto-'jx  integral  by 
the  theory  of  residues. 

It  may  have  been  noticed  that  the  integrals  we  have  evaluated  in  this 
section  have  been  essentially  Fourier  traneforms.  For  example,  (2.?6a)  ahowe  that  the 
Fourier  transform  of  ('k)""^  is  Tr(-i)'^x'^~  sgn  x.  We  shall  verify  this  result  for 
the  cases  n=Oorn=l,  by  using  the  inverse  Fourier  transform  to  show  thrt  the 
Fo-urier  transform  of  sgn  x  is  — j—  and  tlia.t  the  Fourier  transform  of  |i|  is  -  ^  ^ 


nCk)' 


Consider  the  integrals 

I  =  "T-  I        Sgn  X  e   dx 


and 


We  have 


0   2n«i-  00 

-00 


h'ht 


X  e   dx. 
00  '  ' 


I  =  1^  /  sin  kx  dx  =  -^§77-  =  -77- 
0   2tt  «/  ?n'k   TT'k 


and 


0 

.00  -,   J   -00 


It  =  — -  /  X  COS  kx  dx  =  -  -TT—  /   sin  kx  dx  = 5-  , 


which  agrees  with  the  results  above. 

Before  concluding  this  section,  we  wish  to  point  out  thr.t  caution  is 
needed  when  transforming  variables  in  a  finite  part  integral.   Consider 

^  %   =  log  2. 
o 
If  we  make  the  transf  ormp.tion  y  =  ax,  the  Integral  becomes 


/ 
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2a 
r 

0 


7^  =  loF  2a 


instead  of  log  2.   In  order  to  avoid  this  difficulty  we  make  the  following  restriction: 

Restriction:  V/hen  the  varla'hle  oi  integration  in  a  finite  part  inteisral  is  changed . 
the  value  of  the  derivative  of  the  new  variable  with  respect  to  the  ori|&-inal  varia"ble 
must  he  1  at  all  the  singular  points  of  the  integrand. 
As  an  illustration,  consider  the  integral 

/   cot  e  d9  =  log  8inTT/3  =  |  log(3/^). 

Suppose  we  wish  to  make  the  tre,nsformation  u  =  tan  Q/2  in  order  to  obtain  the  inte- 
gral of  a  rational  fijinction.  This  transformation  will  change  the  value  of  the  inte- 
gral. To  satisfy  the  restriction,  we  must  change  the  transformation  to  u=  2tan©/2 
and  then  the  integral  becomes 


/ 


,2(3)"^/^  1  -  ^       .  ,  2 


r^  -^^  =  &og  u  -  i  logd+i-  a 


1  +  SI  ^ 

«  I    log  (3^), 


which  agrees  with  the  original  evalioation. 
Principal  Value   Integrals  and  Singular   Integral  Equations. 
We  have  seen    (page  88)  that 

lim       —-7—    =  -p-  -  iTT6(x). 
a->0    ^^^^  ^ 

Similarly,  we  can  prove  that 


77  =  77  +  i^e(^)- 


lim 

e->0  ^-^^ 

This  result  has   interesting  consequences   in  the   theory  of  analytic  functions.      Suppose 

that  w(z)    is  an  analytic  function  of  z   in  the  upper  half-plane  and  that  w(z)  vanishes 

sufficiently  fast  at   infinity  so  that  an  application  of  Cauchy's   theorem  gives 


1  i°° 


(^(x)dx 
x-z 


-00 

Put  z  =  a+ic  and  consider  the  limit  of  the  integral  as  e  approaches  zero.  We  have 


T  J  1  f      w(y)dr    1  /  w(x)     ^  \      ,y  /  V 

1  ill  -r—r  /     ^  '        =  -r—r  /    , )    '    \dx  +   :r  w  (a )  =  w  (a ) ; 

^„  2ni/    x-a-i€   pTTiy  '(x-a)     2   ^  ' 
€— >0      •'-00 


2ni, 

consequently,  the  value  of  w(z)  on  the  bounctary,  thf-t  is,  the  r'='al  axis,  is  given  by 
the  principal  value  integral 


-Ill  - 

,00 


(2.77)       w(a)=;^/      7^^ 


00 

Now,    suppose  that  u(z)  and  v(z)  are  the  real  and  ima^inpry  parts     respectively,   cf 
the  enalytlt:  fiinction  w(z).      Patting  w  =  u+iv  Into   the  ahove  formula  and  sepsrating 
into  reel  and   iin?ie-lnpry  pprts,  we  have 

(2.78)  ^(«)  =  ny       Tt^.  v(^)  =  -^  y  Tttf^. 

-  00 

These  formulss  show  that  the  resl  part  of  an  analytic  f\inction  can  'he  fo^jnd. 
from  a  knowledge  of  the  imaginary  part^ and  conversely,  the  imaginary  part  can  be 
found  from  the  real  part.  The  methods  of  the  last  section  can  he  used  to  verify  this 
result  in  the  case  of  simple  functions  such  as  w(z)  =  e  . 

The  formulas  (2.78)  can  ^e  written  very  simply  in  terms  of  the  folding  con- 
cept that  was  introduced  on  page  9^.  ''^e  have 

(2.79)  u(x)  =  ^v  *  Cx)"^.  V  =  -  iu  *  Cx)-^  . 

If  we  substitute  the  first  formtda  of  (2.79)  into  the  second  we  get 

V  =  _  ^  (  V  *  Cx)"^)  *  ('x)"^=  -  ^v  •  Cx)"^*  (•x)-^ 

V:  TT 

Since  this  formula  is  valid  for  arbitrary  functions  v(x)  we  must  have 

(2.80)  Cx)"^  *  Cx)"^  =  -  n^  6(x), 

1  2 

or,  expressed  in  words,  the  fold  of  the  principal  value  /x  with  itself  is  -Tr"6(x). 

Another  proof  of  (2.80)  may  be  obtained  by  using  the  Fourier  interral.   By 

Theorem  TV   the  Fourier  transform  of  a  fold  is  the  product  of  the  Fourier  transforms 

of  the  components;  consequently,  since  the  Fourier  transform  of  Cx)   is  ttI  sgn  k 

(see  the  previous  section),  we  find  that  the  Fourier  transform  of  the  fold  of  ( 'x) 

o      2     2 
with  Itself  is  -tt  (sgnk:)  =  -tt  .  This  result  is  obviously  the  Fourier  transform  of 

-TT  P(x),-  therefore  we  again  obtain  (2.80). 

We  shall  now  show  how  (2.80)  may  be  used  to  solve  a  singular  integral 

equation  which  is  important  in  aerodj'namic  theory.   Consider  the  following  integral 

equation  for  the  unknown  fTinction  u(x): 

or 


(2.81)        u(x)  -  \  u*('x)  -^  =  f. 


Take  the  fold  of  this  equation  with  Cx)"  and  use  (2.80).  We  get 
(2.82)        u  *('x)"-'-  +  TT^  \u  =  f  *('x)"-'-. 
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Eliminate  u  *('x)~  between  (2.81)  and  (2.82)  and  we  finally  obtain  the  solution 
of  (2.81): 

f+xf  ♦Cx)"-'- 

¥e  shall  close  this  section  by  deriving  a  formula  which  shows  ths.t  caution 
amst  be  used  when  performing  algebraic  man ipiilat ions  on  principal  value  integrals  in 
more  then  one  dimension.      Consider  the  escpression 


1  fl     ,     l"| 

'(x+y)  L  '^       'yJ   * 


Ordinery  algebraic  manipulations  would   indicate  that  this  expression  equals    ('x)      ('y)~ 
but   this   is  wrong.     Using  the  result    (page  88)   that 

1 


^  €->0      T+€ 


we  have 


(2.83)  =  lim  -rA-y     ^^^^^"f^^     =  li^nl"      ,     f     ,     ,     +  -^  .  -i-1 

^^^^^      (x2+c2)(y2+e2)  [U+e    )(y+c^)       ^^*^         y^+e^J 

again  by  the  results   of  page  88.     T^ormula    (2.8?)   is  useful   in  some  quantum  mechanical 
applications. 

Problems  ,  i    i 

/i""  sinkxdx  /•  °°  e  '^'^' 

2.38  Evaluate  /        i^"*"^     and/  ^     ^  ■  dx. 

-00   ^''"''o^     '^-oo  '(x^-xf) 

2.39  Evaluate     /^  ,,^  J'^ -^^       ifla/bl<l. 

0  0 

2.^0  Make  the  transformation  u  =  2tan(e/2)  in  the  Integrals  of  Problem  2.39  and  then 
evaluate. 


/oo  .  ,  . 
— -  "  „       if  X  has  a  small  positive  imaginary  part.  Oom)a.re  the 


00 

o  x~-x 

0 

sinkxdx 


result  with  the  integral     /       ~     o     o 


0   '(x^-x'^) 

0 


2.^2  Solve  the  Integral  equation 


-  00 
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Tho  discussion  of  the  present  chapter  has  shown  the  utility  and  siorpli- 
city  of  the  theory  of  distrihutions.   It  usually  enables  us  to  obtain  in  a  simple 
way  resulte  which  would  he  fairly  complicated  to  obtain  by  customary  methods.  The 
position  of  the  theory  of  distributions  may  be  compared  to  that  of  the  theory  of 
complex  variables  in  the  stildy  of  acoustic  or  electromagnetic  phenomena.  These 
physical  phenomena  must  be  represented  by  real  numbers,  Howerer,  for  the  purpose 
of  simplifying  the  calculations,  we  introduce  congjlex  numbers  and  at  the  end  inter- 
pret these  numbers  to  obtain  a  real  physical  answer.   In  the  same  way  we  shall 
henceforth,  unless  otherwise  stated,  assume  that  we  are  always  dealing  with  the 
distributions  defined  by  the  functions,  sjnnbollc  or  ordinary,  that  appear  in  our 
calculations.   If  the  final  result  of  the  calculation  is  a  distribution  defined 
"by  an  ordinary  function,  the  ordinary  function  will  be  the  answer  required.  If 
the  final  result  of  the  calculations  Is  a  distribution  defined  by  a  symbolic  function, 
then  the  original  prcblem  does  not  have  a  solution  in  the  usual  sense.  In  this  case, 
however.  If  the  original  problem  Is  re-formulated  so  as  to  apply  to  a  testing  function, 
the  answer  would  be  equal  to  that  obtained  by  applying  the  distribution  defined  by 
the  symbolic  function  to  the  testing  function, 

Wp  shpl''  s^ij!LT.pri7e  somP  of  the  results  of  the  present  c'lriptrr  in  the 
followinf  list  of  •^orTnnlps : 

1.   If  f  (y)  is  p  piecewise  continvouf  f^inction  with  .jur.ps  of  -nrrnitude  a.  ft 

the  T)oint9  p.(  i  =l,?,...v)  and  if  f'C?)  denotes  the  sjmholic  function  for  the 
1 

derivative  of  the  distribution  defined  by  f  (x),  vhile '[f  '  (xlj  is  the  derivative, 

wherever  it  eyists  in  the  ordinary  sense,  of  f (r)j  then 

k 
f(x)  -  CfCxB  +  E  a..  6(x-e^), 


1=1 


EyamTDle.-—  sie-num  x  =  26  (x). 
dx 


2.   If.  fir)   is  a  continuous  diff erentishle  fimction  with  simple  zeros  rt  the 


points   r^,r    .•••»r     then 


n     6(x-rj^) 


HfU))  =  r  if.(r%i  • 

Example.      6(x^-a^)  =  ^- C^  (3c+a)+6  (x-afl 


X     x''6^'^''^x)  =    (-r  ^^^^  e^'^^x).   n  >  0,   m  an  interer. 

n. 

h.    x'^'^e^^^x)  =  0. 

5.      f(x)6(x)  =  f(0)6(x);  g(x)£'(x)  =  -g'(0)6(.v)v(0)£'(x). 


-  llh  - 


iir) 


'■II 


dx  =   lim 
c->0 


y       m 


€         X 


_l-!n 
O.'d-m)     "*■   lJ(2-m) 


i(m-l), 

U-1).' 


7.    /      ^  ^'  '      dr  =   lim 
0    ('x)P  ^->0_ 


IT.   integral 

0       .P 


/        _p     ^-^        O.'(l-p)      *•••*       (m-l)i 


m-p 


p  non-intepral;   m  is   the   lerf-eat   integer  contained   in  p. 
8.     ^^^  =  t.^^=-t-i-6(^). 


dx 


dx 


(m), 


9.  f-  C('xr^H(72  =  -inCxr'^-^HCxHc-r  ~-^ . 


^  [:('^rPH(xa  =  -p('xr^-^F(x). 


dx 
d 


dx: 


10 


11.     u  = 


.Hkl 


(n-2)('x)"-2 


is  p.  solntion  of  (x   -a')'  =  -6. 


12.      lim 


2      2 
c->0  e  +x 


n 


=  tt6(x).    lim-^=  -^ 
e  +x 

n^ (n) I 


k  e       dk  =   2tt(.-i;    u        Ix;, 


-00 
.00 


,v^  .,  .n+1    , 

Ik,      f       k^e^^dk  =  TT(-i)^6(:,K  ^^ ^,  . 

>/o  Cx)^-*"^ 

15.  /       sfD  k  e       dk  =  ~  . 
-co 

16.  6(x)6(y;  =  r;^  -^        r     =  x  +:•   . 

6(r-r   )6(0) 

17.  f(x)6(y)e(.-7    )  =  r . 


7y^T 


6(x)e(y)f(^) 


N   _iki 


i.  Tvr 


2  • 


2  2 .    2^   2  ,    ft      o 

r     =   /  +v  +z,7=rsinO,    2     =r. 
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